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SECTION  I 


INTRODUCTION 

The  typical  Auxiliary  Power  Unit  (APU)  is  a  small  gas  turbine,  independent 
of  the  main  engine,  which  supplies  shaft  horsepower  and  compressed  air  to  start 
the  main  engine  and  run  accessories  as  needed  in  the  aircraft.  APU  designs  are 
typically  single-spool  or  single-spool  with  power  turbine.  The  need  for  high 
power  density  in  a  confined  space  necessitates  use  of  both  axial  and  centrifugal 
designs  in  the  compressors  and  turbines.  Presently  there  exists  a  need  to 
analyze  APU's  and  aerodynamic  torque  converter  configurations  in  AFAPL. 

This  report  serves  to  document  the  previous  modifications  to  the  SMOTE 
program  which  have  resulted  in  the  APUCODE  program.  APUCODE  simplifies  the 
SMOTE  program  to  analyze  the  basic  single-spool  gas  generator  while  incorporating 
options  peculiar  to  APU's.  At  present,  APUCODE  analyses  only  design-point 
conditions.  In  addition,  numerous  changes  have  been  made  to  make  input  and 
output  pages  more  readable. 

This  report  is  a  user's  manual  for  design-point  APUCODE.  It  includes 
complete  Fortran  listing  of  the  code  and  sample  data  cases.  The  code  has  been 
programmed  in  a  straightforward  manner  for  simplicity  and  ease  of  modification. 
Comment  cards  describe  the  major  subdivisions  of  the  main  subroutine,  APU,  and 
flowcharts  aid  the  understanding  of  the  program  logic. 
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SECTION  II 


PROGRAM  DESCRIPTION 


1 .  APU  IN  GENERAL 

APU  is  the  main  program  subroutine  in  APUCODE.  It  controls  the  use  of 
subroutines,  sequences  the  logic  of  the  calculations,  and  prints  the  output. 

The  program  logic  follows  the  gas  path  through  the  APU  from  inlet  to  exhaust. 
Thermodynamic  properties  of  the  flow  are  computed  at  the  inlet  or  exit  of  each 
component.  Typical  computing  stations  are  shown  in  Figure  1  and  defined  in 
Table  1.  Calculations  are  assumed  to  be  adiabatic.  The  properties  at  down¬ 
stream  computing  stations  are  calculated  from  mass-averaged  upstream  quantities, 
and  the  change  in  those  quantities  due  to  addition  or  extraction  of  shaft  work 
and  flow. 

2.  SEQUENCE  OF  OPERATIONS  IN  APU 

a.  Inlet  and  Compressors 

After  printing  the  INPUT  page,  inlet  conditions  are  computed  from  1962 
standard  atmosphere.  Inlet  duct  pressure  loss  is  used  to  compute  the  total 
pressure  at  the  compressor  inlet.  Throughout  the  program,  at  each  computing 
station,  subroutine  THERMO  provides  calculation  of  gas  properties.  Work  for  the 
low  and  high  pressure  compressors  is  computed  from  input  pressure  ratios,  effi¬ 
ciencies,  and  losses.  If  only  one  compressor  exists,  it  will  be  computed  as  the 
High  Pressure  Compressor  (HPC). 

b.  Bleed  Flows  from  the  High  Pressure  Compressor 

Bleed  flows  for  turbine  cooling  and  accessory  power  are  subtracted  from 
the  HPC  inlet  flow.  The  HPC  bleed  flow  can  be  split  into  turbine  inlet  nozzle 
cooling.  High  Pressure  Turbine  (HPT)  cooling,  and  Low  Pressure  Turbine  (LPT) 
cooling.  The  remainder  is  bleed  air  which  is  lost  from  the  APU  flow  and  used 
to  power  the  APU  fuel  pump,  oil  pump  and  other  accessories.  Each  bleed  flow  is 
accounted  for  as  a  fraction  of  the  HPC  inlet  flow  and  listed  in  the  program  output. 
Losses  in  the  bleed  ducting  are  used  to  arrive  at  the  thermodynamic  properties 
of  the  bleed  flow  when  it  enters  the  main  flow  to  serve  as  turbine  cooling  air. 

c.  Combustor 

The  combustor  calculates  the  properties  at  combustor  exit  of  a  mixture 
of  air  and  standard  JP-4  fuel.  The  combustor  exit  temperature,  T4,  must  be  input. 
Fuel  flow  and  fuel  heating  value  are  optional  and  will  be  calculated  if  not 
specified.  The  combustor  calculation  Iterates  to  a  final  exit  temperature  using 


TABLE  1 


COMPUTING  STATION  DESCRIPTION 


ambient  conditions 

low  pressure  compressor  inlet 

low  pressure  compressor  discharge 

high  pressure  compressor  inlet 

high  pressure  compressor  discharge/burner  inlet 

main  combustor  discharge 

high  pressure  turbine  inlet  plus  inlet  nozzle 
cooling 

high  pressure  turbine  discharge  before  cooling 
flow  has  been  added 

low  pressure  turbine  inlet  plus  HPT  cooling 

low  pressure  turbine  discharge  before  cooling 
flow  has  been  added 

low  pressure  turbine  discharge  including  LPT  cooling 
exhaust  duct  exit 


the  Air  Force  Quadratic  Interpolation  Routine  (AFQUIR).  If  the  calculation  fails 
to  converge,  a  message  is  printed  to  that  effect, 
d.  Turbines 

There  are  three  choices  for  the  turbine  configuration,  see  Figure  2:  1)  a 

one-stage  turbine,  2)  a  two-stage  turbine,  or  3)  a  one-stage  turbine  followed  by 
a  free  turbine.  The  program  makes  two  passes  during  each  case:  the  first  employs 
the  pressure  ratios  and  efficiencies  as  input  to  calculate  the  shaft  work  extracted 
from  the  APU;  the  second  pass  calculates  the  turbine  pressure  ratios  necessary  to 
provide  the  design  APU  shaft  horsepower.  At  the  end  of  the  first  pass,  the  output 
page  is  printed  and  control  is  transferred  back  to  the  combustor  calculation. 

With  compressors  unchanged,  the  program  continues  on  the  second  pass,  and  prints 
the  output  for  the  second  pass. 

(1)  Turbine  Overexpansion 

If,  on  the  first  pass,  the  last  turbine  lowers  the  turbine  exit  pressure 
to  less  than  atmospheric,  an  error  is  printed  indicating  overexpansion  and  the 
turbine  exit  pressure  is  set  equal  to  atmospheric.  Then,  on  the  second  pass, 
pressure  ratio  in  the  last  turbine  is  recomputed  based  on  expansion  to  atmospheric 
turbine  exit  pressure. 

(2)  Turbine  Cooling 

Before  output  is  printed,  cooling  air  modifies  the  first  turbine 
inlet  nozzles  and  each  turbine  exit.  It  is  assumed  that  the  cooling  air  injected 
into  the  turbine  flow  does  no  useful  work,  thus  cooling  air  is  added  to  the  flow 
computation  after  the  appropriate  turbine  exit  station.  It  is  then  assumed  to 
mix  completely  with  the  main  flow  and  thus  perform  useful  work  in  the  downstream 
turbine. 

(3)  Exit  Velocity 

Flow  velocity  is  computed  at  the  last  turbine  exit  based  on  the 
change  in  enthalpy  from  turbine  exit  to  ambient  conditions.  Thrust  and  SFC  are 
computed  from  this  velocity.  On  the  first  pass,  when  turbine  overexpansion  occurs, 
this  enthalpy  change  will  be  negative;  hence,  the  exit  velocity  calculation  will 
be  invalid.  This  velocity  and  all  values  dependent  on  it  will  be  printed  out  as 
"I." 

(4)  Peculiarities  in  the  Case  of  One  Turbine  Followed  by  a  Power  Turbine 

Because  the  power  turbine  rotates  a  shaft  separate  from  the  compressor- 

turbine  shaft,  the  power  turbine  requires  special  attention.  The  power  derived 
from  the  main  turbine  must  be  exactly  equal  to  the  power  required  by  the  compressors 


2  STAGE  COMPRESSOR  1  STAGE  TURBINE  WITH  POWER  TURBINE 
CONFIGURATION 


Figure  2  -  Possible  APU  Configurations 


and  the  accessories.  Thus,  on  the  first  pass,  if  the  input  main  turbine  pressure 
ratio  results  in  this  requirement  not  being  satisfied,  the  correct  pressure  ratio 
of  the  main  turbine  is  calculated.  For  this  reason,  the  first  pass  output 
pressure  ratio  for  the  main  turbine  may  not  be  that  which  was  input. 

Continuing  in  the  first  pass,  the  power  turbine  work  is  calculated 
from  the  input  pressure  ratio  and  printed  out.  In  the  second  pass,  the  main 
turbine  is  unchanged,  and  the  power  turbine  pressure  ratio  is  recalculated  to 
provide  the  design  shaft  horsepower. 

3.  QUADRATIC  INTERPOLATION  ROUTINE 

Throughout  the  program  there  are  many  small  loops  (for  example,  thermodynamic 
iterations  and  table  lookup)  which  require  convergence.  Trial -and-error  methods 
and  linear  interpolations  can  be  time-consuming,  especially  when  a  tight  tolerance 
is  necessary;  therefore,  a  general  interpolation  routine  called  AFQUIR  was 
developed. 

a.  Variables 

This  routine  requires  a  dummy  array  dimensioned  for  nine  locations.  The 
other  input  variables  are  the  independent  and  the  dependent  variables,  the  answer 
or  value  which  the  dependent  variable  is  to  converge  upon,  the  number  of  tries 
at  convergence,  the  tolerance,  and  a  variable  called  DIR. 

b.  The  DIR 

The  DIR  is  either  set  or  calculated  in  the  calling  program.  It  is  an 
initial  guess  at  the  direction  and  percentage  change  applied  to  the  first  value 
of  the  independent  variable.  If  not  enough  is  known  to  calculate  a  DIR,  an 
arbitrary  value  may  be  set.  This  should  not  affect  the  final  result  but  may 
increase  the  number  of  tries  at  convergence. 

c.  The  Logic 

The  DIR,  thus,  establishes  the  second  value  of  the  independent  variable. 
This  value  is  used  in  the  calling  program  to  determine  a  corresponding  second 
value  of  the  dependent  variable.  AFQUIR  is  called  a  second  time  with  these  two 
sets  of  values.  A  linear  interpolation  results  in  a  third  value  of  the  independent 
variable.  AFQUIR  is  then  called  a  third  time  with  all  three  pairs  of  independent 
and  dependent  variables  and  a  quadratic  interpolation  is  made.  The  values  of  these 
three  pairs  of  variables  are  stored  in  the  dummy  array  and,  henceforth,  quadratic 
interpolations  are  made  using  the  three  sets  which  give  values  closest  to  the 
answer.  Values  farthest  from  the  answer  are  lost. 
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d.  Safeguards 

Various  safeguards  are  built  into  AFQUIR  which  either  return  the 
interpolation  method  to  the  initial  guess  or  linear  interpolation  if  the  roots 
of  the  quadratic  become  complex.  Complex  roots  are  encountered  if  the  quadratic 
does  not  intercept  the  answer,  if  the  value  of  the  independent  variable  differs 
radically  from  previous  values,  or  if  two  pairs  of  independent  and  dependent 
variables  are  identical.  It  is  also  possible  to  preload  the  dummy  array  and  to 
start  directly  at  the  linear  or  quadratic  interpolations  if  desired. 

e.  Summary 

In  summary,  AFQUIR  is  a  completely  flexible  routine  which  performs 
quadratic  interpolation  for  quick  convergence  of  general  functions.  (SMOTE  manual, 
reference  no.  1,  pages  11  and  12). 

4.  THERMODYNAMIC  PROPERTIES  SUBROUTINES 

a.  THERMO  and  PROCOM 

Two  subroutines  are  used  to  calculate  thermodynamic  properties  of  the 
flow  at  each  computing  station.  They  are  PROCOM,  which  calculates  gas  properties 
for  either  air  or  fuel-air  mixture,  and  THERMO,  which  calls  PROCOM  and  calculates 
more  complex  properties.  These  subroutines  are  fairly  simple  and  are  listed  with 
the  rest  of  the  code  in  Appendix  B. 

b.  Use 

The  use  of  PROCOM  is  straightforward,  but  THERMO  may  be  used  in  several 
ways.  If  pressure  and  enthalpy  are  input,  total  temperature  is  calculated.  If 
pressure  and  temperature  are  supplied,  enthalpy  is  calculated.  Fuel-to-air  ratio 
may  also  be  input  and  taken  into  account  in  the  property  computations.  These 
options  are  set  by  the  programmer  and  are  not  available  as  options  to  the  user 
from  the  input  deck. 


SECTION  III 


INPUT  DESCRIPTION 


1 .  GENERAL 

Subroutine  APU  is  the  main  subroutine  of  APUCODE.  Its  logic  is  straightforward 
and  is  described  at  major  subdivisions  by  comment  cards.  At  the  beginning  of 
the  code  listing,  the  coirment  cards  list  brief  descriptions  of  the  program 
variables.  The  complete  set  of  input  variables  need  not  be  input  for  any  one 
case.  Variables  not  input  are  set  equal  to  their  default  values  as  listed  in 
Section  1 1 1. 6.  Those  parameters  which  are  input  override  the  default  values. 

2.  CARD  FORMAT 

All  of  the  data  cards  must  start  in  Column  2  or  beyond.  To  start  a  data  case, 
punch  $INPUT  in  Columns  2-7  then  skip  a  space  and  continue  by  punching  the  input 
parameter  names,  followed  by  an  "equals"  sign  and  the  appropriate  numerical  value. 
Values  can  be  input  up  to  column  72  inclusively.  There  is  no  limit  to  the  number 
of  cards  in  a  data  case.  Finish  the  data  case  by  punching  a  $-sign  immediately 
following  the  last  value  in  that  data  case.  The  program  executes  the  data  case 
and  returns  to  read  the  next  case.  The  program  will  use  values  from  the  preceding 
point  unless  specifically  changed  in  the  data  case  that  follows.  See  example  A, 
Setup  of  APUCODE  Deck,  in  Appendix  A.  To  terminate  a  string  of  data  cases,  the 
last  data  case  should  be:  $INPUT  I END  =  1$. 

3.  TITLE  CARD 

If  specified,  a  title  will  be  printed  at  the  top  of  each  case.  To  read  a 
title  card,  set  ITITLE  =  1  in  the  NAMELIST  data  to  which  the  title  applies.  The 
program  then  reads  Columns  1-66  of  the  first  card  following  that  data  case  as  the 
title.  Since  the  program  automatically  sets  ITITLE  =  0  after  a  title  has  been 
read,  ITITLE  has  to  be  set  =  1  only  when  a  new  title  is  desired. 

4.  ESTIMATES  FOR  ONE  CASE 

To  compile  (CP  time)  *  7.2  sec. 

Input/Output  time  (10)  =  23  sec. 

Line  Estimate  =  255  lines  per  case,  or  3  computer  pages 

5.  SUBROUTINE  DESCRIPTION  AND  PROGRAMMER  OPTIONS 

a.  ATMOS  (ALTP,  T1 ,  PI)  -  Computes  temperature  and  pressure  of  the  1962 
U.  S.  standard  atmosphere  up  to  90  kilometers. 

b.  PR0C0M  (FAR,  T1 ,  CS,  AKEX,  CP,  R,  PHI,  H)  -  Calculates  thermodynamic 
gas  properties  for  either  air,  or  a  fuel -air  mixture  based  on  JP-4. 


c.  THERMO  (£,  H,  T,  SOUT,  AM,  CP,  6A,  L,  FAR,  K)  -  Calculates  thermodynamic 
conditions  based  on  PROCOM.  If  K=l,  P,  FAR,  and  H  are  input;  if  K=0,  T,  FAR, 
and  P  are  input. 

d.  COMB  (PIN,  T4,  HIN,  ETAB,  DPB,  POUT,  HOUT,  SOUT,  FAR,  HV)  -  Computes 
stoichiometric  temperature  in  combustor  if  heating  value  for  fuel  and  fuel  flow 
are  not  input. 

e.  C0MB2  (PIN,  T4,  HIN.,  ETAB,  DPB,  POUT,  HOUT,  SOUT,  FAR,  HV)  -  Similar  to 
COMB.  Used  when  fuel  heating  value  is  not  input  but  fuel  flow  is  input. 

f.  AFQUIR  (X,  AIND,  DEPEND,  ANS,  AJ,  TOL,  DIR»  ANEW,  ICON)  -  General  quadratic 
interpolation  routine.  All  parameters  explained  in  computer  listing,  see  AppendixB. 
NOTE:  Parameters  which  are  input  to  the  subroutine  are  underlined. 

6.  INPUT  DESCRIPTION 

The  following  list  describes  all  possible  input  parameters  and,  in  parenthesis, 
their  default  values  if  they  are  omitted: 

ALTP  -  ALTITUDE  IN  FEET  (0.0) 

DELTAMB  -  CHANGE  IN  AMBIENT  TEMPERATURE  FROM  1962  ATMOSPHERE  (0.0) 

DELPAMB  -  CHANGE  IN  AMBIENT  PRESSURE  FROM  1962  PRESSURE  (0.0) 

DPBLEED  -  BLEED  DUCT  PRESSURE  LOSS  (0.0) 

DPB  -  MAIN  COMBUSTOR  PRESSURE  LOSS  (P3-P4)/P3  (0.0) 

DPC  -  INTERCOMPRESSOR  PRESSURE  LOSS  (P22L-P22)/P22L  (0.0) 

DPD  -  INLET  DUCT  PRESSURE  LOSS  (0.0) 

DPT  -  INTERTURBINE  PRESSURE  LOSS  (P5-P51 )/P5  (0.0) 

ETALPC  -  LP  COMPRESSOR  EFFICIENCY  (1.0) 

ETAHPC  -  HP  COMPRESSOR  EFFICIENCY  (1.0) 

ETAB  -  MAIN  COMBUSTION  EFFICIENCY  (1.0) 

ETAHPT  -  HP  TURBINE  EFFICIENCY  (1.0) 

ELALPT  -  LP  TURBINE  EFFICIENCY  (1 .0) 

ETAM  -  MECHANICAL  EFFICIENCY  (1.0) 

*  HV  -  MAIN  COMBUSTOR  FUEL  HEATING  VALUE  AT  T4  FOR  JP-4 

(Optional.  If  omitted,  it  is  calculated  by  the  combustor 
calculation.) 

HPDS  -  DESIGN  SHAFT  POWER  (0.0) 

HCSL  -  HPC  SEAL  LEAKAGE/HPC  FLOW  (0.0) 

I END  -  0  A  DATA  CASE  FOLLOWS 

1  THIS  IS  THE  LAST  CASE  (0) 
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ITITLE  -  0  NO  TITLE  CARD  WILL  BE  INPUT 

1  TITLE  CARD  FOLLOWS  THE  DATA  CASE  (0) 

LCSL  -  LPC  SEAL  LEAKAGE/LPC  FLOW  (0.0) 

M  -  NUMBER  OF  TURBINE  STAGES  (no  default) 

=  1,  ONE  TURBINE  APU  DESCRIBED  BY  HPC  INPUT 
=  2,  TWO  TURBINE  APU 

=  4,  ONE  TURBINE  DESCRIBED  BY  THE  HPT  INPUT  AND  ONE 
POWER  TURBINE  DESCRIBED  BY  LPT  INPUT 
N  -  NUMBER  OF  COMPRESSOR  STAGES  (no  default) 

=  1,  ONE  COMPRESSOR  STAGE  DESCRIBED  BY  HPC  INPUT 
=  2,  TWO  COMPRESSOR  STAGES 

PRES  -  A  MULTIPLIER  ON  PRESSURE  VALUES  FOR  OUTPUT  (1.0) 

=  14.7,  PRESSURES  ARE  IN  PSI 
=  1.0,  PRESSURES  ARE  IN  ATMOSPHERES 

*  PRHPT  -  HP  TURBINE  PRESSURE  RATIO  (1.0) 

*  PRLPT  -  LP  TURBINE  PRESSURE  RATIO  (1.0) 

PCBLT4  -  HP  TURBINE  INLET  NOZZLE  COOLING  FLOW/HP  COMPRESSOR  INLET  FLOW  (0.0) 
PCBLHP  -  HP  TURBINE  EXIT  COOLING  FLOW/HP  COMPRESSOR  INLET  FLOW  (0.0) 

PCBLLP  -  LP  TURBINE  EXIT  COOLING  FLOW/HP  COMPRESSOR  INLET  FLOW  (0.0) 

PRLPC  -  LP  COMPRESSOR  PRESSURE  RATIO  (1.0) 

PRHPC  -  HP  COMPRESSOR  PRESSURE  RATIO  (1.0) 

*  T4  -  COMBUSTOR  DISCHARGE  TEMPERATURE  (DEGREE  RANKINE)  (must  be  input, 

no  default) 

WAENG  -  TOTAL  ENGINE  INLET  AIRFLOW  (1.0) 

WF4DS  -  DESIGN  FUEL  FLOW  (LB/HR)  (if  omitted.it  will  be  calculated  in 
combustor  calculations) 

The  program  recalculates  those  parameters  preceded  by  an  asterisk.  For  this  reason, 
they  must  be  input  for  each  case  if  the  user  desires  to  keep  them  the  same  for  each 


SECTION  IV 


OUTPUT  DESCRIPTION 


1.  REGULAR  OUTPUT 

Output  for  the  three  typical  APU  configurations  are  presented  in  Example 
Cases  A,  B,  and  C  of  the  Appendix  A.  All  descriptors  on  the  printed  output 
are  designed  to  be  self-explanatory. 

a.  Two  Compressor  -  One  Turbine  APU 

Example  A  shows  an  APU  with  two  compressor  stages  and  one  turbine 
stage.  Fuel  flow  and  heating  value  are  input.  Pass  1  output  shows  that 
using  the  input  parameters,  the  shaft  horsepower  obtained  exceeds  the  desired 
design  shaft  horsepower.  Pass  2  recalculates  the  pressure  ratio  of  the  turbine 
to  satisfy  the  design  shaft  horsepower. 

b.  Two  Turbine  APU 

If  two  turbine  stages  had  been  used,  the  program  would  have  recomputed 
both  HPT  and  LPT  pressure  ratios  to  obtain  the  design  horsepower  by  preserving 
the  work-split  in  the  turbines  as  it  was  in  Pass  1. 

c.  Example  B 

Example  B  illustrates  the  calculation  of  fuel  heating  value  when  fuel 
flow  is  input. 

d.  Power  Turbine  APU 

Example  C  exhibits  a  power  turbine  configuration  with  fuel  heating 
value  Input  and  fuel  flow  calculated.  The  input  HPT  pressure  ratio  was  deficient 
so  HPT  pressure  ratio  was  recalculated  on  Pass  1.  With  the  input  power  turbine 
pressure  ratio,  overexpansion  occurred.  Exit  velocity  and  its  dependent  param¬ 
eters  were  computed  as  being  indefinite  errors  and  printed  as  "I."  Pass  2 
remedied  the  problem  by  recalculating  power  turbine  pressure  ratio  to  achieve 
design  shaft  horsepower. 

2.  DIAGNOSTICS 

Input  errors  will  be  indicated  for  each  component  when  parameters  are  not 
logical.  For  compressors  and  turbines,  execution  stops  when  an  input  pressure 
ratio  is  less  than  1.0  or  an  efficiency  is  greater  than  1.0.  For  the  combustor, 
execution  stops  when  input  combustor  exit  temperature  or  fuel  flow  Is  equal 
to  or  less  than  zero,  or  when  combustor  efficiency  is  greater  than  1.0. 
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a.  PROCOM  Diagnostics 

Subroutine  PROCOM  generates  diagnostics  In  four  cases:  1)  when  fuel- 
air  ratio  Is  above  0.067623,  2)  fuel-air  ratio  Is  below  zero,  3)  the  temperature 
input  to  PROCOM  Is  below  300°R,  or  4)  the  temperature  Input  to  PROCOM  is  above 
4000° R. 

b.  THERMO  Diagnostics 

Subroutine  THERMO  halts  execution  after  15  Iterations,  If  the  value  of 
enthalpy  Is  not  within  a  tolerance  of  0.001  percent  of  the  previous  enthalpy 
guess. 

c.  Combustor  Diagnostics 

Subroutine  COMB  and  C0MB2  both  print  a  message  indicating  lack  of 
convergence  of  stoichiometric  temperature.  This  message  Is  a  direct  result  of 
failure  in  subroutine  AFQUIR.  If  after  15  Iterations  the  solution  to  the 
quadratic  equation  for  stoichiometric  temperature  has  not  been  found,  execution 
stops. 
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APPENDIX  A 
EXAMPLE  CASES 


INPUT  CARDS  FOR  EXAMPLE  CASE  A 


IINPUT  ALTP=0. ,  DELTAMB=71 . ,  PRES=14.696,  WAENG-2. 485346,  ETAB-.995, 
DPC=.01 ,  PRLPC=1 .794,  PRHPC=6.0,  LCSL*0.0,  HCSL=.02,  ETALPC«.78,  DPB*.07, 
ETALPT=.870,  T4=2559.67,  DPT*0.0,  M«2,  N*2,  WF4DS=0. ,  ETAHPC*.746, 
HPACC=9.5,  PRLPT=9.418,  HPDS®200. ,  HV=18400. ,  ITITLE=1 ,  DPIN=.04. 
ETAM=.98,  ETAHPT=.8,  PRHPT=1.5,  DPBLEED* . 01 4$ 

SAMPLE  DATA  FOR  APU 
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SAMPLE  DATA  FOR  APU 


♦  ♦♦INPUT  VALUE  S ♦♦ ♦ 


ALTITUDE  (FT)  *  0. 

PRESSURE  (PSI)  *  14.696 

DELTA  TEMP  FROM  1962  ATMOS  -  71.00000 

BEGIN  INDICATOR  -  1 
END  INDICATOR  -  0 

NO  OF  COMPRESSOR  STAGES  -  2 
NO  OF  TURBINE  STAGES  «  l 

FLOW  *  2.48636 

MFCH  EFF  *  .98000 

HP  ACC  EXTRACTED  FOR  ACCESS  *  9.60 

OESIGN  SHAFT  POWER  *  166.0000 


INLET  DUCT  PRESS 

LOSS 

.04 

LPC 

PRESS  RATIO 

-  1.794 

LPC 

EF.  =.  .  78000 

LPC 

SEAL  LEAKAGE 

*  0. 

0 

LPC 

BLEED »  TO  BE 

USEO 

LATER 

INTER-COMP  PRESS 

LOSS 

M 

.01 

BLEED  DUCT  PRESS 

LOSS 

m 

.014 

HP  C 

PRESS  RATIC 

»  6.000 

HP  C  EFF  *  .74600 

HP C  SFAL  LFAKAGE  *  .02 

MAIN  COMB  PRESS  LOSS  *  .07000 

MAIN  COMBUST  I  ON  EFF  =  .99600 

OESIGN  FUEL  FLOW  (LB/HR)  -  160.00 

MAIN  COM3  FUEL  HEATING  VALUE  AT  T4  FOR  JP4  «  18400. 
COMB  DISCHARGE  TEMP  (DEGREE  RANKlNE)  -  2669.67 
HPT  INLET  NOZ  COOLING/HPC  INLET  FLOW  -  0.00000 

HPT  PRESS  RATIO  *  1.000 
HPT  EFF  *  .80 

TURB  COOLING/HPC  INLET  FLO*  «  0.00000 

INTER  TURB  PRESS  LOSS  (P5-P8D/P5  -  0.00 

LPT  PRESS  RATIO  -  9.418 
LPT  EFF  «  .87000 

EXHAUST  DUCT  PRESS  LOSS  (PSJ-P1)/PS3  -  0.0 
LPT  COOLING/HPC  INLET  FLOW  *  .01000 


Example  A  -  2  Compressor  Stages,  1  Turbine  Stage 
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♦♦♦INPUT  VALUES*** 


ALTITUDE  TFT)  *  0. 

PRESSURE  (PSl)  -  LA. 696 

DELTA  TEMP  FROM  1962  ATMOS  -  71.00000 

BEGIN  INDICATOR  *  0 
ENO  INDICATOR  *  0 

NO  OF  COMPRESSOR  STAGES  *  2 
NO  OF  TURBINE  STAGES  «  l 

FLO*  *  2.A8535 

MECH  EFF  *  .98000 

HP  ACC  EXTRACTED  FOR  ACCESS  *  9.50 

DESIGN  SHAFT  POWER  »  165.0000 

INLET  DUCT  PRESS  LOSS  -  .OA 
LPC  PRESS  RATIO  *  1.79A 
LPC  EFF  =  .78000 

LPC  SEAL  LEAKAGE  *  0.0 


LPC 

BLEED, 

TO  BE 

USED 

LATER 

INTER-COMP 

PRESS 

LOSS 

-  .01 

BLE 

ED  DUCT 

PRESS 

LOSS 

*  .  01 A 

HPC 

PRESS 

RATIO 

»  6.000 

HP C  EFF  =  • 7A600 

HPC  SEAL  LEAKAGE  -  .02 

MAIN  COMB  PRESS  LOSS  -  .07000 

MAIN  COMBUSTION  EFF  *  .99500 

DESIGN  FUEL  FLO*  (LB/HR)  *  190.00 

MAIN  COM3  FUEL  HEATING  VALUE  AT  TA  FOR  JPA  *  0. 

COMB  DISCHARGE  TEMP  (DEGREE  RANKINE)  »  2559.67 
HPT  INLET  NOZ  COOLING/HPC  INLET  FLOW  -  0.00000 

HPT  PRESS  RATIO  *  1.000 
HPT  EFF  -  .80 

TURB  COOLING/HPC  INLET  FLOW  ■  0.00000 

INTER  TURB  PRESS  LOSS  (P5-P5U/P5  -  0.00 

LPT  PRESS  RATIO  »  9.A18 
LPT  EFF  *  .87000 

EXHAUST  DUCT  PRESS  LOSS'  (P53-PD/P53  *  0.0 
LPT  COOLING/HPC  INLET  FLOW  *  .01000 


Example  B  -  2  Stage  Compressor,  1  Stage  Turbine 
With  Fuel  Heating  Value  Calculated 
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EXPENDABLE  CAS  IF  IE k  JF S 


*** INPUT  VALUE  S  *  *  * 


ALTITUDE  (FT)  =  0. 

PRESSURE  (PSD  =  14.696 

DELTA  TtMP  FROM  1962  ATMOS  -  0.00000 

BEGIN  INDICATOR  =*  1 
END  INDICATOR  «  0 

NO  OF  COMPRESSOR  STAGES  =  1 
NO  OF  TURBINE  STAGES  =  4 

FLCW  *'  3.80200 

MECH  EFF  =  . 980C0 

HP. A CC  EXTRACTED  FOR  ACCESS  *  4.00 

DESIGN  SHAFT  POWER  =  230.0000 

INLET  DUCT  PRESS  LOSS  =  .01 

LPC  PRESS  RATIO  =  1.000 
LPC  EFF  =  I.OCOOO 
LPC  SEAL  LEAKAGE  =  0.0 
LPC  BLEED i  TO  BE  USED  LATER 

INTER-COM?  PRESS  LOSS  =  0.00 
3LFED  DUCT  PRESS  LOSS  =  0.000 

HPC  PRESS  RATIO  «  2.827 
HPC  EFF  =  .79200 

HPC  SFAL  LEAKAGE  =  0.00 

MAIN  COMB  PRESS  LOSS  =  .12300 

MAIN  COMBUSTION  EFF  =  .95000 

DESIGN  FUEL  FLOW  (LB/HR)  =  0.00 

MAIN  COMB  FUEL  HEATING  VALUE  AT  T4  FOR  JP 4  =  18400. 
COMB  DISCHARGE  TEMP  ( DEGREE  RANKINE)  =  2260.00 
HPT  INLET  NO 7  COQLING/HPC  INLET  FLOW  *  O.OOCOO 

HPT  PRESS  RATIC  =  1.540 
HPT  FFF  =  .83 

TUR3  COOLING/HPC  INLET  FLOW  «  0.00000 

INTER  TURB  PRESS  LOSS  (P5-P51I/P5  *  .02 

LPT  PRESS  RATIO  =  I  .600 
LPT  EFF  *  .85000 

EXHAUST  DUCT  PRESS  LOSS  (P53-PD/P53  -  0.0 
LPT  COOLING/HPC  INLET  FLOW  *  O.COOOO 


Example  C  -  1  Stage  Compressor,  1  Stage  Turbine  with  Power  Turbine 


♦♦♦OUTPUT  VALUIS*** 
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PROGRAM  APU  (INPUT,  OUTPJT,  TAPE5  =TNPUT , 

$  T AP£6=0UT°UT) 

COMMON  /P  ARAM/  ETA. PC,  ETAH^C,  ETAB,  ETAHPT , 

$  FTALPT,  PRLPC, 

1  PRHPC,  DPB,  DPT,  Q?E,  0=>IN,  D®C,  HPACC,  H^DS, 

$  LCSL ,  HCSL , 

2  PCBLC ,  T4,  ETAM 

DIMENSION  ANUM  (2),  AWORO  (2),  3WORO  (2),  TITLE 
$  (11) 

REAL  LCSL 

DATA  AWORD  /  3HATM,  3H:,SI  / 

DATA  ANUM  /  3H  ,  3H  2  / 

DATA  BWORD  /  3HLP  ,  6HHPACC  / 

C 

NAMELIST  /  INPUT  /  ALTP,  DE.PAMB,  DELTAMB,  PRES, 

5  WAENG,  PRHPT, 

1  0P3LEE0,  DPIN,  OPG,  PR.3C,  ^PH^C,  LCSL,  HCSL, 

$  ETALPC,  ETAHPC, 

2  ETAHPT,  LTALPT ,  T4,  DPT ,  M,  HV,  N,  ETAB,  WF4DS, 

$  DPB,  PCULT4, 

3  PC BL HP,  FTAM,  PCBLLP,  ITITLE,  »’LPT,  HPDS,  HPAC3, 

$  IEND 

C 

G  STATION  1  -  AMBIENT  CQNDI1IONS 

G  STATION  2  -  LP  COMPRESSOR  INLET 

G  STATION  21  -  LP  COMPRESSOR  DISCHARGE 

G  STATION  22  -  HP  COMPRESSOR  INLET 

C  STATION  3  -  HP  COMPRESSOR  DISCH ARGE/PURNtR  INLET 

G  STATION  4  -  MAIN  COMBUSTOR  DISCHARGE 

C  STATION  41  -  HP  TURBINE  INLET 

C  STATION  5  -  HP  TURBINE  DISCHARGE  BEFORE  COOLING  FLOW 

C  $  HAS  BEbN  ADDED 

C  STATION  51  -  LP  TURBINE  INLET 

C  STATION  52  -  LP  TURBINE  OISCHAPGF  BEFORE  COOLING  FLOW 

G  $  HAS  BEEN  AOOEO 

C  STATION  53  -  LP  TURBINE  OISGHARGE  INCLUDING  COOLING 

C  $  FLOW  EFFECTS 

G  STATION  6  -  EXHAUST  DUCT  EXIT 

C  ALTP  -  ALTITUDE  IN  FEET 
G  DPBLEF.D  -  BLEED  DUCT  PRESSURE  LOSS 

C  BLTOTL  -  HP  COMPRESSOR  TOTAL  BLEED  (FRACTION  OF  WA22) 
C  BLT4  -  HP  TURBINE  INLET  NOZZLE  COOLING  FLOW  (FRACT  OF 
C  $  BLTOTL) 

C  BLHPT  -  HP  TUR8INE  OISGHARGE  COOLING  FLOW  (FRACT  OF 
C  $  BLTOTL) 

C  BLLPT  -  LP  TURBINE  DFSCHAQGE  CDOLING  FLOW  (FRACT  OF 
C  $  BLTOTL) 

C  CS  -  AMBIENT  SPEFD  OF  SOUND  <rPS) 

C  OELTAMB  -  CHANGE  IN  AMBIENT  TEM»ERATJRF  FROM  1962 
G  S  ATMOSPHERE 
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o  n  n  o  o  o  u  u  u  o  u  o  o  u  o  o  u  u  n  a  o  u  u  o  o  o  o  (i  o  o  u  o  o  o  o  u  u  u  o  o  o  o  u  o  u  o  u  o  u  o 


DZLPAMB 

0P9 

OPC 

DPIN 

OPE 

OPT 

FAR 

ETALPC 

ETAHPC 

tTAB 

ETAHPT 

ETALPT 

ETAM 

ETATHM 

FN/WA 

FG 

FN 

HS6 

HV 

$ 

H  PL  PC 
HPHPC 
HPACC 
HPM 
$ 

H»MAMB 

$ 

HPtXT 
HPDS 
HCSL  - 
LCSL 
M 


$ 

N 

POT 

PRES 

PROA 

PRHPT 

PRLPT 

PC9LA 

S 

PCBLT4 

$ 

PC9LHP 

$ 


-  CHANGE  IN  AMBIENT  PRESSURE  FROM  1962  PRESSURE 
MAIN  COMBUSTOR  PRESSURE  .OSS  (P3-P4)/P3 
INTER-COMPRESSOR  PRESSURE  LOSS  ( P22L-P22) /P22L 

-  INI  FT  DUCT  PRESSURE  LOSS 

EXHAUST  DUCT  PRESSURE  LOSS  (P53-PD/P53 
INTER-TURBINE  PRESSURE  .OSS  (P5-P5U/P5 
FUEL  FLOW/AIRFLOW 

-  LP  COMPRESSOR  EFFICIENCY 

-  HP  COMPRESSOR  EFFICIENCY 
*AIN  COMBUSTION  EFFICIENCY 

-  HP  TURBINE  EFFICIENCY 

-  LP  TURBINt  EFFICIENCY 
MECHANICAL  EFFICIENCY 
THERMAL  efficiency 

NET  THRUST/TOTAL  ENSINr  AIR-.OM 
GROSS  THRUST 
NET  THRUST 

STATIC  ENTHALPY  AT  EXHAJST  EXI"  ( BTU/LB  MASS) 
MAIN  COMBUSTOR  “UEL  HEATING  YALUS  AT  T4  FOR  JP- 
4 

HORSEPOWER  NEEDED  FO?  LP  CON»RFSSOR 
HORSEPOWER  NEEDED  POR  HP  COMPRESSOR 
-  HOPStPOWER  EXTRACTED  FOR  ASSESSORltS 
SHP  EXTRACTED,  CORRECTED  FOR  MECHANICAL 
EFFICIENCY 

SHP  EXTRACTED,  CORRECTED  FOR  S.L.  ,  STANDARD  DAY 
CONDITIONS 

SHAFT  HORSEPOWER  EXTRACTED 
DFSIGN  SHAFT  POWER 
HPC  SEAL  LFAKAGE'H°C  -10* 

LPC  SEAL  LEAKAGE/LPC  "LOW 
NUMBER  OF  TURBINE  STAGES 
=  i,  ONc.  TURBINE  APJ 
=  2,  TWO  TUR3INE  APJ 

=  4,  TWO  TURBINE  4®U ,  ONE  Or  WHICH  IS  A  POWER 
TURBINE 

NUMBER  OF  COMPRESSOR  STAGES 

OESIGN  POINT  WORK  LO&D  ’-RCENTAGE  PER  TURBINE 
A  MULTIPLIER  ON  PRESSURE  VA.UES  FOR  OUTPUT 
OVEPALL  CORE  COMPRESSION  RATIO  (P3/P2) 

-  HP  TURBINE  PRESSURE  RATIO 

-  LP  TURBINE  PRESSURE  RATTO 

ACCESSORY  AIR  B*£ED  rLOW  EXTRACTED  FROM  HP 
COMPRESSOR/ 

HPC  INLET  FLOW 

HP  TURBINE  INLET  NOZZLE  COOLING  FLOW  /  HP 
COMPRESSOR 
INLET  FLOW 

HP  TURBINt  EXIT  COOLING  rLOW  /  HP  COMPRESSOR 
INLET  FLOW 
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c 


c 

c 

c 

L 

r \ 
-J 

C 

c 

c 

-/ 

c 

c 

c 

c 

c 

w 

c 


PC3LLP 

prJpc 

PP.HPC 

SFCU 

$ 

T4 

$ 

TS6 

THG6 

V6 

WA22 

WAENG 

WF 

WF4DS 

WG 

WLPC 

WHPC 


-  LP  TUkBINt  EXIT  CODLING  “LOW  /  H°  COMPRESSOR 

INLET  FLOW 

-  LP  COMPRESSOR  PRESSURE  RATIO 
-  HP  COMPRESSOR  PRESSURE  PATIO 

-  UNINSTALLED  SPECIFIC  FUEL  CONSUMPTION  (LB  FUEL/ 

HR- HP) 

-  COMBUSTOR  DISCHARGE  TEM3ERATJRE  (DEGREE 

RANKINF) 

-  STATIC  TEMPERATURE  AT  EXHAUST  EXIT  (DEGRtE  F) 

-  COkE  GROSS  THRUST  PER  PDUNO  TOTAL  AIRFLOW 

-  EXHAUST  EXIT  VE-OCITT 

-  AIRFLOW  AT  HPC  INtET/WAENG 

-  TOTAL  ENGINE  INLET  AIRF.OW 

-  FUEL  FLOW  RATE  (LB/SEC) 

-  DESIGN  FUEL  FLOW  UB/HR) 

-  MASS  FLOW  RATE  (INCLJDES  FUEL  FLOW)  (LB/SEC) 

-  WORK  OF  THE  LP  COMPPESSDR  (3TU/LB) 

-  WORK  OF  THE  HP  COMPRESSOR  (BTIJ/LB) 


***  CONSTANTS 

DATA  G,  AJ,  I END j  IP  /  32.174)49,  778.26,  0 ,  1  / 
PRES  =  WAFNG  =  PRL3C  =  PRHPC  =  ETALPC  =  ETAHPC 

$  =  1. 

ETAHPT  =  FTALPT  =  ETAM  =  £TAT  =  °RHPT  =  PRLPT  = 


$ 

1. 

WF4DS 

=  OELTAMR  =  D£L3AM3  =  CPIN  *  DPB 

PPC 

=  DPT 

=  DPI  =  .:SL  =  HC5L  =  PCBLC 

$ 

=  0, 

PCBLLP 

=  PCBLHP  =  HPACC  =  HPDS  =  OPBLEFD 

s 

0. 

IPAS 

=  0 

IN 

=  IB  = 

2 

READ  (5, 

INPUT) 

IF  ( I END 

•  ED.  1) 

WRITE  (6,  1)0) 

IF  (IEND 

•  ECU  1) 

STOP 

IF  (ITITLE  -  1) 

2 ,  3 ,  3 

WRITE  (6, 

101) 

GO  TO  4 

READ  (5, 

102  ) 

(TITLE(I),  1=1,  11) 

WRITE  (6, 

103) 

(TITLE ( I)  ,  1=1,  11) 

WRITE  (6, 

104) 

WRITE  (6, 

10  5) 

ALT3 

WPITE  (6, 

106) 

PRES 

WRITE  (6, 

107) 

DELTAMB 

WRITE  (6, 

108) 

ITITLE 

WRITE  (6, 

109) 

IEND 

WPITE  (6, 

110) 

N 

WRITE  (6, 

111) 

M 

WRITE  (6, 

112) 

WAENG 

0. 

PCBLT4 
HV  = 
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o  o 


HRI  r  E 

(6, 

113) 

ETAM 

WRITE 

<G» 

114) 

HPADC 

WRITE 

(fa, 

115) 

HP03 

WRITE 

(b, 

116) 

DPH 

WRITE 

<6, 

117) 

PRL3C 

WRITE 

<6, 

118) 

etalp: 

WRITE 

<6> 

119) 

LCSL 

WPIT  £ 

(6, 

120) 

WRITE 

<6, 

121) 

OPC 

WRITE 

(6, 

122) 

DPB_FED 

WRITE 

<6, 

123) 

PRH3C 

WRITE 

(6, 

124) 

E  T  AH°C 

WRITE 

(6, 

125) 

HCS. 

WPITE 

(6* 

126) 

0P9 

WHITE 

(fa, 

127) 

ET  A3 

WRITE 

(6, 

128) 

WF4DS 

WRITE 

<6, 

129) 

HV 

WRITE 

(6, 

130) 

T  4 

WRITE 

<6, 

131) 

PCBlT  4 

WRITE 

<6, 

132) 

PRHaT 

WRITE 

<6, 

133) 

ETAHPT 

WRITE 

(6, 

134) 

PCB-HP 

WRITE 

(6, 

135) 

DPT 

WRITE 

(6, 

136) 

PRL’T 

WRITE 

(6, 

137) 

ETAlPT 

WDITE 

(b, 

138) 

OPE 

WPITE 

(6, 

13  9) 

PC9.LP 

WRITE 

(6, 

140) 

CALL  ATMOS 

(A  LTP 

,  Til,  P 

Ti  =  Til  +  DELTAMB 

Pi  =  Pi  +  DELSAMB 

CAcL  PROCOM  (0.0,  Ti,  CS,  XI,  X2,  XT,  X4,  X5) 

C 

C  BEGINNING  OF  INLET 

CALL  THERMO  (Pi,  Hi,  TI,  31,  AMHl,  C°1 ,  GA1 ,  3,  O.P, 
$  0) 

H2  =  Hi 

P2  =  Pi  *  (i.a  -  Q  =  IN> 

CALL  THERMO  (  P2 ,  H2,  T2,  S2,  AMW2,  CP2 ,  GA2,  0,  J.P, 
$  1) 

HPM  =  0.0 

IF  ( PRLPC  .FQ.  1.0)  30  TO  34 

IF  ( PP.L  PC  .LT.  1.0  .OR.  ETAL3D  .5T.  1.0)  GO  TO  35 
IF  ( PRHPC  .LT.  1.0  .DR.  ETAHPC  .GT.  1.0)  GO  TO  35 
5  IF  (N  .EQ.  1)  GO  TO  3 

***  BEGINNING  OF  LP  COMPRESSOR 
P21  =  PRLPC  *  3 2 

PATIO  A  =  <  GA2  -  1.)  /  3 A  2 

WL PC  =  ( CP2  *  T?  »  <<(P2i  /  ?2)  **  RATIOA)  - 

$  1.))  /  ETALPO 
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OO  OO 


H21  =  WLPC  «■  H2 

CAl.  THERMO  (P2i,  H21,  T 21,  S21,  AMW21,  CP21,  GA21, 

%  u  »  0.0,  1) 

GO  TO  7 

6  P21  =  P2 

T2i  =  T  2 

H21  =  H  2 

WLPC  =  0.0 

7  IF  (PRHPC  .t0.  1.0)  GO  TO  8 

***  BEGINNING  OF  HP  COMPRESSOR 

P22  =  P 21  *  (1.0  -  5PC) 

H22  =  H  21 

T  22  =  T  21 

CALL  THERMO  (P22,  H22  ,  T22,  S?2,  AMW22 »  ^22,  &A22, 
$  0,  0.0,  0) 


P3 

= 

PRHPC 

*  322 

KATIOB 

= 

(GA22 

-  1. )  /  GA?? 

WHPC 

= 

(  CP22 

*  T 22  *  U(P3  /  P?2>  **  RATIOS)  - 

$ 

1.))  /  tTAHPC 

H3 

= 

WHPC  *■ 

H2? 

CALL  THERMO 

<P3,  H3 ,  T 3 ,  S3,  AMW3 ,  CP3 ,  GA 3 ,  3, 

J.  0, 

s 

1) 

RATI03 

s 

(GA3  - 

l.J)  /  G A3 

GO  TO  9 

8 

P22 

= 

P  21 

P3 

= 

P  22 

T  3 

= 

T  21 

H  3 

= 

H  21 

WHPC 

i).C 

CALL  THERMO 

( P3 ,  H3 ,  T3,  S3,  AMW3 ,  CP3,  GA 3 ,  0, 

0.3, 

$ 

0) 

9 

PROA 

= 

PRHPC 

*  PRLPC  *  (1.0  -  n°C) 

•  * 

AIR  FLOWS, 

FRACTIONS  OF  WAING, 

WA2 

= 

l.C 

WA21 

= 

WA2  * 

(1.0  -  _CS.) 

WA22 

= 

WA21 

WAIRH® 

= 

CP3  * 

T3  *  ( ( ( 33  t  PI)  **  RATIOS)  - 

1.  0) 

BLA 

= 

(HPAlC 

*  .  706  B>  f  WAIRH® 

PCBLA 

s 

BLA  / 

(WA22  *  WA*MG> 

PCBLC 

s 

PCBLT4 

♦  PCOLWP  P03LLP  ♦  PCBLA 

WA3 

s 

WA22  * 

(1.0  -  HC3-  -  PCBLC) 

IF  ( ( PCBLC  .LT.  1,0))  GO  TO  13 

PRINT  *,  ••PERCENTS  O"  BLEED  rL0W3  .GT.  1.0  — HALT 
J  EXECUTION" 

GO  TO  1 

10  BLT4  =  PCBLT4  *  WA22 
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BLHPT  =  PC8LHP  *  W422 

6LLPT  *  PCBLLP  *  M A22 

BLP  s  P3  *  (1.3  -  O’BLIED) 

HPLPC  s  (WLPC  *  ^ A 2  *  WAEMG)  /  .7368 

HPHPC  =  ( WHPC  *  M  A22  *  WAENG)  /  .7068 

TB  =  T  3 

CALL  THFPHO  (BLP,  H8,  TB,  SB,  AMHB,  CRB,  GAB,  0, 

S  L  .  0  ,  0  ) 

***  BEGINNING  OF  MAIN  COMBUSTOR 

IF  (WF40S  .LT.  0.0  .OR.  T4  ,.T.  3.0  .OR. 

$  E.TAB.  GT.l.  0  )  GO  TO  33 

11  IF  (WF4DS  .GT.  (5.0)  GO  TO  14 
IF  (HV  .GT.  0.0)  GO  TO  12 

CALL  COMB  <P3,  T4,  H3,  FT AB ,  DP3,  R4,  H4,  S4,  FAR4, 

$  HV) 

GO  TO  13 

12  P4  =  P3  *  (1.0  -  03B  > 

FAR4  =  (.276  *  (T4  -  T3»)  t  ( (ETA8  *  HV)  - 

$  (.27b  *  ( T4  - 
1  T3) )  ) 

CALL  THERMO  ( P4 ,  H4,  T4,  34,  AMW4,  CP4 ,  GA4,  1, 

I  FAR4,  G> 

FAR4  =  ( CP4  *  <T4  -  T3>)  /  ( (ETAB  *  HV)  -  (CP4 

$  *  (T4  -  T 3 ) )) 

13  WF4  =  FAR4  *  WA3 

GO  TO  15 

14  WF4  *  V4F4DS  /  (W4ENG  *  3600.) 

F AR4  *  WF4  /  WA3 

P4  =  P3  *  (1.)  -  03B) 

CALL  THERMO  (P4,  H4,  S4,  AMW4 »  CP4 ,  GA4,  1, 

$  FAR4 ,  0) 

IF  (HV  .GT.  0.)  GO  TO  15 

IF  (HV  .FQ.  O.u)  HV  =  (  (C3*  *  (T4  -  T3)  /  FA34>  ♦ 

$  CP4  *  (T4  - 
1  T3 ) )  /  ETAB 

CALL  C0M82  (P3,  T4,  H3,  ETAE,  DP3,  *4,  H4,  S4,  FAP4 , 
$  HV) 

15  WG4  =  W A3  ♦  HF4 

P41  =  P4 

WG41  =  WG4  +  BLT4 

FAR41  *  HF4  /  (WA3  RLT4) 

H41  =  ( H4  *  WG4  *  HB  *  BLT4)  /  WG41 

CAlL  THERMO  (P41,  HhI,  T41,  S41,  A MR41 ,  CP41,  GA41, 

S  1,  FAR41,  1) 

IF  ( 8LT 4  ,LE.  0.3)  T41  =  T4 

IF  ( PRHPT  .LT.  1,0  .JR.  ETAH3T  • GT.  1.0)  GO  TO  37 
IF  (PKLPT  .LT.  1.0  .OR.  ETALPT  .GT.  1.0)  GO  TO  37 
IF  (M  .tQ,  1)  GO  TO  16 
IF  (M  •  EQ.  2)  GO  TO  19 
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IF  (M  • EQ.  4)  GO  TO  21 

***  ONE  TURBINE  APU 
16  P5  =  P41 

PRHPT  =  P41  /  PS 

H5  =  H  41 

T5  =  T  41 

FAR5  =  FAR41 

CALL  THtPMO  (P5»  H5>  T5»  55,  AMW5,  CP*>,  GA5» 

$  FAR?  ,  1) 

P61  =  P  5 

HE1  =  H5 

T51  =  T5 

WG51  =  WG41 

FAR51  =  FAR41 

CALL  THFPMO  (P51»  Htil ,  T51,  S51,  AMW51,  0P51 ,  GA51 , 

$  1,  FAR51,  1) 

IF  (HPM  .NE.  0  •  U )  GO  TD  17 
RATIO  =  (GA51  “  1  •  0 >  *  GA51 

TWORK  =  ETALPT  *  0*51  *  T51  *  (1.0  -  ((1.0  / 

$  PRLPT)  ** 

1  RATIO)) 

DHTL  =  TWORK 

H52  =  H51  -  DHTL 

P52  =  P51  /  PRi.PT 

HPLT  =  (TWORK  *  WG51  *  HAEN3)  /  .7*68 

HP EXT  =  HPLT  -  H^C  -  HPHPC 

HPM  =  (HPtXT  -  HPACO  »  ETA* 

PCT1  =  0.0 

PCT2  =  1.0 

IPAS  =  1 

GO  TO  23 

17  HPBTU  =  (  (HPACC  (HP 3S  f  E.TAN)  )  *  .7C68)  / 

%  ( WG51  *  WAENG) 

OHTL  =  HPBTU  ♦  ( ( <H3  -  H’?>  *  WA22  +  (H21  -  H2) 

$  *  WA2)  / 

1  WG51 ) 

HB2  =  H51  -  DHTL 

P52  =  051  *  (1.0  -  (DHT.  f  (ETALPT  *  HD51  * 

%  T 51)  ) )  ♦*  (GA51 
1  /  (GA51  -  1.0)) 

CALL  THERMO  (P52,  H52,  T52,  352,  AMWi>2,  CP52,  GA52 , 
$  1,  FAR51,  1) 

TWOKK  =  OHTL 

MP.T  =  (TWORK  *  WS51  *  <AEN3)  /  .7168 

HPEXT  =  HPLT  -  HPLPC  -  HPHPC 

PCT1  =  0.0 

PCT2  =  1.0 

HPM  =  (HPtXT  -  H*AC:)  *  FTAM 

IPAS  =  2 
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***  TWO  TURBINE  APU 

Id  CALL  THERMO  (PL  1,  H41,  TM  ,  541,  AMW41,  C^4l,  6A41, 
S  1,  FAR41,  0) 

IF  ( HPM  .NE.  0.0)  GU  T 3  13 
RATIO  =  (GA41  -  1.0)  /  GA41 

THWORK  =  ETAHFT  *  GP41  *  T  41  *  (l.J  -  (C1.0  / 

$  PRHPT)  ** 

1  RATIO)) 

DHTH  =  THWORK 

H5  =  H41  -  DHTH 

P5  =  P41  ✓  PRHPT 

CALL  THERMO  <P5,  H5,  T5,  55,  A«M3,  Ct>5,  GA5,  1, 

$  FAR.41,  1) 

P51  =  P5  *  (1.0  -  D*T) 

WG51  =  WG41  +  dLHPT 

H51  =  (H5  *  WG41  t  (B  »  BL HPT)  /  WG51 

FAR51  =  WF4  /  ( W351  -  WF4 ) 

CALL  THERMO  (P51,  H51 ,  T51,  S51,  AMW51 ,  CP51,  GA51 , 
$  lt  FAR51,  1) 

RATIO  =  ( GA51  -  1.0)  /  GA51 

DHTL  =  ETALPT  *  C»51  *  T51  *  (l.i  -  ((1.0  / 

$  PRLPT)  ** 

1  RATIO)) 

H52  =  H 51  -  DHTL 

P52  =  P51  /  PRLPT 

CALL  THERMO  (P52,  H52,  T52,  552,  AMW52,  CP52,  5A52, 
$  I*  FARF1,  1) 

PCT 1  =  DHTH  /  (DHTH  f  DHTL) 

PCT2  =  DHTL  /  (DHTH  «•  DHTL) 

IPAS  =  1 

HPHT  =  (THWORK  *  W341  *  WAEMG)  /  .7068 

HPLT  =  (DHTL  *  HG51  *  WAEMG)  /  .7068 

HPEXT  =  HPHT  -  H=»HPC  H»„T  -  H°LPC 

HPM  =  (HPEXT  -  H®AC  3)  *  ETAM 

GO  TO  23 

19  HPPTU  =  ((HPACC  *>  (HP3S  /  ETAM)  )  *  .7168)  / 

$  (WG41  *  WAENG) 

DHALL  =  HPBTU  ♦  (  ( H3  -  H22)  *  WA22  ♦  (H21  -  H2) 

S  *  WA2)  /  WG41 

DWORKHT  =  DHALL  *  (P3T1)  ♦  WG41  *  WAENG 

FAR5  =  FAR41 

DHTH  =  DWORKHT  /  (WGM  *  WAIMG) 

H5  =  H41  -  DHTH 

P5  =  P41  *  (1,0  -  (DHTH  /  (FTAHPT  *  0P41  * 

$  T 41)  )  )  **  (GA41 
1  i  (GA41  -  1.0)) 

CALL  THERMO  (  P5 ,  H5 ,  T5,  55,  AMW5 »  OP5 ,  GA5,  1, 

$  FAR5,  1) 
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PRHPT  =  P41  /  P 5 

Pi51  =  P5  *  (1.0  -  Q3T) 

WG51  =  WG41  «•  B'.HPT 

H51  =  (H5  *  WG41  «•  HB  *  RLHPT)  /  WGEi 

FAR51  =  WF4  /  ( W35 1  -  WF4 ) 

CALL  THLRMO  (P51,  H3l,  Tbl,  551,  AMW51,  CP51  ,  GA51, 

$  1,  FAR51,  1) 

DWORKLT  =  DHALL  *  (PRLPT  /  ( PRHPT  *  PRLPT) )  *  WG51 
$  *  WAENG 

DHTL  =  DWORKLT  /  (WG51  *  WAENG) 

H52  =  H 51  -  DHTL 

P52  =  Pbl  *  (l.u  -  OHT.  /  (ETALPT  *  OP51  * 

3  T51)  ) )  **  (GA51 
1  /  (GA51  -  1.0)) 

CALL  THEPMO  ( P52 ,  H52,  T52,  S52,  AMW52,  CP52,  GA52, 

$  1,  FAR51,  1) 

PCT 1  =  DWORKHT  /  (DWORKHT  ♦  DWORKLT) 

PCT2  =  DWORKLT  /  (DWDRKHT  +•  DWORKLT) 

HPHT  =  DWORKHT  /  .7038 

HPLT  =  DWORKLT  /  .7058 

HPEXT  =  HPHT  -  HPHDC  HPLT  -  HPLPO 

HPM  =  ( HPFXT  -  HPACC)  *  ET AH 

IPAS  =  2 

GO  TO  23 

***  TWO  TURBINE  APU  WITH  ONE  0“  THE  TURBINES  A  POWER 
$  TURBINE 

20  CALL  THERMO  (P41,  H4l,  TH,  $41,  AMW41,  CP41  ,  GA41, 

$  1,  FAR41,  b ) 

RATIO  =  ( GA41  -  1.0)  f  GA41 

THWORK  =  ETAHPT  *  CP41  *  T41  *  (1.0  -  ((1.0  / 

$  PRHPT)  *♦ 

1  RATIO)) 

OHTH  =  THWORK 

H5  =  H41  -  OHTH 

P5  =  P4l  /  PRHPT 

CmLL  THERMO  ( P5 ,  H5,  T5,  S5 ,  AMW5 ,  CP5,  GA5,  1, 

$  FAR41,  1) 

HPHT  =  (THWORK  *  WG41  *  H AEHG)  /  .PDE8 

HPTT  =  HPHT  -  H3ACD  -  HPHP3  -  HPLPO 

IF  (HPTT  .GE.  0.0)  GO  TO  ?1 

HP3TU  =  ( HP ACC  *  .7053)  7  (W341  *  WAENG) 

OHTH  =  HP8TU  +  (((H3  -  H22)  *  WA22  ♦  (H21  -  HP) 

S  *  WA2)  / 

1  WG41) 

H5  =  H 41  -  DHTH 

P5  =  P41  *  (1.0  -  (DHTH  /  (ETAH»T  *  CP41  * 

i  T41) )  )  *•  ( GA 41 
1  f  (GA41  -  1.0) 

CALL  THERMO  (PS,  H5,  T5,  55,  AHW5 ,  CP5 ,  GAS,  1, 

S  F AR4 1 ,  1) 


48 


o  o  o 


PRHPT  =  P41  /  P5 

PRINT  *,  "PP.E SSURE  RATIO  CALO JLATEO  AS  NEEDED  FOR 
$  HPT" 

21  PR1  =  P5  •  (1.3  -  0»T) 

WG51  a  WG41  ♦  3,-HPT 

H51  =  (H5  *  WG41  ♦  HB  *  BLHPT  >  /  WGS1 

FAR51  =  WF4  /  (W351  -  WF!») 

CALL  THERMO  ( P51 ,  H51,  T51,  S51,  AMM51,  CP51,  GA51, 

$  1,  FARSI,  1) 

IF  ( HPM  »NL.  O.H)  GO  TO  22 
RATI02  =  (GA51  -1.0)  *  GA51 

TWORK  =  ETALPT  *  C»51  *  T51  *  (1.0  -  ((1.0  / 

$  PRLPT)  ♦* 

1  RATI02) ) 

OHTL  =  TWORK 

H52  =  H 51  -  OHTL 

PS2  =  P51  /  PRLPT 

CALL  THERMO  (P52,  H52,  T52,  S52,  AMM52,  CP52,  GA52, 

$  1,  FAR51,  1) 

HPEXT  =  (TWORK  ♦  HG51  *  W  A  FN  3 )  /  .7068 

HPLT  a  HPEXT 

HPM  a  HPEXT  *  £T AM 

OWORKHT  =  OHTH  *  W341  *  WAENG 

DWORKLT  a  OHTL  *  W351  *  WAENG 

PCT 1  *  OWORKHT  /  ( DWORKMT  +  DWORKLT) 

PCT2  *  DWORKLT  /  (OWORKHT  +  DWORKLT) 

IPAS  a  l 

GO  TO  23 

22  OHTL  a  ( (HPDS  /  ETAMJ  *  .7058)  /  (WG51  *  WAENG) 

CALL  THERMO  (P52,  H52,  T52,  S52,  AMW52,  <^52,  GA52, 

S  1»  FAR51,  1) 

P52  a  P51  *  (1.0  -  (OHTL  f  (ETALPT  *  CP51  * 

$  T51) ) )  **  ( GA51 
1  /  (GA51  -  1.0)) 

H52  a  H51  -  DHTL 

TWORK  a  DHTL 

HPEXT  a  (TWORK  *  W351  *  HAEN3)  /  .7068 

HPLT  a  HPEXT 

OWORKHT  a  OHTH  *  WS41  *  WAENG 

DWORKLT  a  OHTL  *  WG51  *  WAENG 

PCT1  a  OWORKHT  /  (OWORKHT  *  DWORKLT) 

PCT2  =  DWORKLT  7  (DWORKHT  ♦  DWORKLT) 

HPM  a  HPEXT  *  ETAM 

PRLPT  a  P51  /  P5? 

IPAS  a  2 

GO  TO  23 

**•  LOW  PRESSURE  TURBINE  OR  POWER  TUR3INE  EXIT,  WITH 
$  COOLING  Flow, 
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IF 

ANY,  PLUS  STATION  5  CALCULATTQN. 

23 

PRLPT 

=  P51  /  P52 

P52 

=  P 51  /  PR. PT 

CALL  THERMO  (P52,  H52,  M2,  552,  AMW52, 

CP52 ,  GA52, 

S 

1,  FAR51,  1) 

P53 

=  P  52 

IF  (IPAS  . EQ.  1)  GO  TO  2 7 

IF  (P53 

•GE.  Pi)  GO  TO  2 7 

ON 

2ND  PASS  WHEN  WORK  REQUIRED  CAUSES 

$ 

OVEREXPANSION  OF  .6  ST 

TURBINE,  ASSUME  THE  TJRBIM5”  FXPANOS  ONLY  TO 

$ 

ATMOS pHERIC. 

PRINT  * 

,  •*P53<Pi,LPT  OVEREXPANDED,  SET 

P53=P1, 

$ 

RECAL3 

WORK" 

P53 

=  P  52  =  Pi 

reprtl 

=  P5i  /  PS? 

PTTIO 

=  ( GA51  -  i.)  /  GA51 

REWORK 

=  ETALPT  *  CB51  *  T51  *  (1.  - 

((1.  / 

$ 

REPRTL) 

**  RFTIO)) 

PRLPT 

=  REPRTL 

H52 

=  H5i  -  REWORK 

HPRE 

=  (REWORK  *  WG5 1  *  WAEN3 )  f  • 

7068 

IF  (M  . 

GT.  i)  GO  TO  24 

HPEXT 

=  HPRE  -  HbLP3  -  HPHPC 

HPM 

=  (HPEXT  -  HPAC3)  *  EMM 

GO  TO  26 

24 

IF  (M  , 

GT.  2)  GO  TO  25 

HPEXT 

=  HPHT  -  HPHPC  f  H»RE  -  H=»LPC 

HPM 

=  (HPEXT  -  HPACC)  *  EMM 

PCT2 

=  H PRC  /  (HPRE  *■  H>HT) 

GO  TO  26 

25 

HPEXT 

=  HPRE 

PCT2 

=  HPRE  /  (HPHT  «■  HP ?E) 

HPLT 

=  HPEXT 

HPM 

=  HPEXT  ♦  ETAM 

26 

CALL  THERMO  (P52,  H52,  TB2,  552,  AMW52 , 

CP52,  GA57, 

$ 

1,  FAR51,  1) 

ON  1ST  PASS,  IF  LAST  TURBINF  NOT 

$ 

OVEREXPANDEO, 

27 

BS6 

=  Pi 

LAST  TURBINE  EXIT  COOLING  CALCULATIONS 

WG53 

=  WG51  «■  81LPT 

H53 

=  (H52  *  W351  HB  *  3LLPT)  / 

W653 

FAR53 

=  WF4  /  (W353  -  WFM 

GALL  THFRMO  (P53,  H53,  T53,  5*53,  AMW53, 

CP53,  GA53, 

$ 

1,  FAR53,  1) 

RTI053 

=  ( GA53  -  1.)  t  GA5 3 

OPE 

=  ( P53  -  P56)  t  °53 

TS6 

=  T  53 

CALL  THfRMO  (PS6,  HS3 ,  TSb,  SS6,  AMW6,  Cp6 ,  GAb,  1, 
S  PAR53,  U) 
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c  *** 


c 

c  *** 


DELNOZ 

r 

H53  *  (1.  -  <31  f  P53)  **  RTI053) 

V6 

s 

SQRT (2  *  G  *  (DELNO?)  *  AJ) 

WG6 

= 

WG53 

FAR5 

= 

FAR53 

Pfc RFORMANCE 

CALCULATIONS 

THG6 

V 6  *  MG6  /  G 

F  NOW  A 

= 

THG6  /  WA2 

SFCU 

r 

WF4  *  WAENG  *  3633  /  HPM 

FG 

= 

THG6  *  WAENG 

FN 

s 

FNQWA  *  WAENG 

ETATHM 

s 

HPlXT  *  S50.  /  ( Wr4  *  WAENG  *  HV 

OLFINE  OUTPUT  PARAMETERS 
XWA21  =  WA21  *  WAENG 

XWA22  =  WA22  *  WAENG 

XWA3  =  WA3  *  WAENG 

X  WG4  =  WG4  *  WAENG 

XWG41  =  WG41  *  WAENG 

XWG51  =  WGi>i  *  WAENG 

XHG53  =  WG53  *  WAENG 

X WG6  =  WG6  *  WAENG 

XWF4  =  WF4  *  WAINS  *  3635. 

X°1  =  PI  *  PRES 

XP2  =  P2  *  PRES 

XP21  =  P21  *  PRES 

XP22  =  P 22  *  PRES 

XP3  =  P3  *  PRES 

XP4  =  P4  *  PRES 


XP41 

XP5 

XP51 

XP52 

XP53 

XPS6 

CF2 

PI) 

CF21 

CF22 

CF3 

CF4 

CF41 

CF51 

CF52 

CF53 

CFo 

CFF 

HPMAMB 

THETA1 


P41  *  PRES 
P5  *  PRES 
P 51  *  PRES 
P52  *  PRES 
P53  *  PRES 
PS6  *  PRES 

( WA2  *  WAENG  *  S1RT(T2  /  Til))  /  C  P2  / 

(XWA21  *  SXRTIT21  /  Til ) )  /  (P21  /  PI) 

(XWA22  *  SORT  <  T  22  t  Til))  /  CP2?  /  PI) 

(XWA3  *  S0RTU3  f  Til))  /  <P3  /  PI) 

CXWG4  *  SORT ( T4  f  Til))  /  IP4  t  PI) 
(XWG41  *  SORT ( T 4L  f  Til))  /  <P41  /  PI) 
(XWG51  *  SORT  <  T51  f  Til))  /  (P51  /  PI) 

( XWG51  *  SORT  € T 52  /  Til))  /  CP52  /  PI) 

( XWG53  *  SORT (T53  f  Til))  /  (P53  /  PI) 
CXWG6  *  SQRTUS6  f  Til))  7  (PS6  /  PI) 

X WF4  /  <<P4  /  PI)  *  SORT (T4  /  Til)) 

< HPEXT  *  SQRTCTi  f  518.67))  /  (»1  /  1.0) 
T 1  /  Til 
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THETA2 

r 

T2  /  Til 

THETA3 

s 

T  3  /  Til 

THtT  A21 

= 

T 21  /  Tit 

THETA22 

9 

T 22  /  Til 

THETA4 

= 

T4  /  Til 

THETA41 

9 

T 41  /  Til 

THETA5 

9 

T 5  /  Til 

THETA51 

= 

T 51  /  Til 

TH£T  Ac2 

X 

T 52  /  Til 

THETA53 

s 

T  F1  /  Til 

THETA6 

r 

'  S6  /  Til 

OELTA1 

s 

*1  /  PI 

0ELTA2 

= 

P2  /  PI 

OFLTA21 

- 

P21  /  PI 

OELTA22 

= 

P  22  /  PI 

DELTA3 

s 

P 3  /  »i 

0LLTA4 

s 

P4  /  PI 

DELTA41 

s 

P41  /  PI 

DELTAS 

9 

P 5  /  PI 

DELTA51 

r 

P51  /  PI 

DLLTA52 

= 

P52  /  PI 

DELTA53 

9 

P53  /  PI 

Dfc  LTA6 

= 

PS6  /  PI 

Tl 

9 

Tl  -  459.87 

T2 

= 

T 2  -  459.67 

T  21 

= 

T  21  -  459.57 

T22 

s 

T  22  -  459.57 

TB 

- 

TB  -  459.57 

T  3 

9 

T  3  -  459.67 

T4 

= 

T 4  -  459.67 

T  41 

= 

T 41  -  459.57 

T5 

9 

T 5  -  459.67 

T51 

= 

T  51  -  4*3.57 

T52 

= 

T 52  -  459.57 

T53 

X 

T 53  -  459.57 

TS6 

r 

TS6  -  453,67 

***  WRITE  OUT  THE  OUTPUT  3AGF 

IF  (PRES  .GT.  1.0)  IP  =  2 
IF  ( ITITLc  «cQ.  w)  GO  TO  28 
WRITE  (6,  141)  AWORO  <  T  3) 

GO  TO  29 

28  WRITE  <6,  142)  AWORD  (I3) 

29  ITITLE  =  0 

WRITE  (6,  143)  XP1,  Tl,  HI,  Di.TAl,  THETAt 

IF  (PRLPC  .LE.  1.0)  30  TO  30 

WRITE  (6,  144)  XP2,  T2,  H2,  WAENG,  CF2,  0ELTA2, 

%  THETA2,  ETALPC, 

1  3RLPC,  DPIN 

WRITE  (6,  145)  XP21,  T21,  -121 ,  XWA21,  CF21, 

$  0ELTA21,  THETA21, 
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1  DPC 

3T  IF  (PRHPC  .LE.  l.u)  10  TO  31 

WRITE  (b,  146)  XP22,  T22,  422,  XWA22,  CF22, 

S  0LLTA22,  THETA22, 

1  ETA4PC,  PRHPC 
31  IF  (PRHPC  .LE.  1,0  13  =  1 

WRITE  (6,  147)  flWOlO  (13),  XP 3,  T3,  H3,  XWA3 ,  Cr3, 

$  0ELTA3, 

1  THETA3,  PROA 

WRITE  (6,  148)  XP4,  T4,  H4,  XW34,  CF4,  DELTA4, 

$  THETA4,  FAR4, 

1  ETAQ,  OPB 

IF  ( H  .EQ.  1)  GO  TO  32 

WRITE  (6,  149)  XP41,  T41,  441,  XWG41,  CF41, 

$  DLLTA41,  THETA41, 

1  FAR41,  ETAHPT,  PRH3T 

WRITE  (6,  150)  XP5,  T5,  H5,  XW341,  CF41,  DELTAS, 

$  THETA5,  FAR41  , 

1  OPT 


WRITE 

(6, 

151) 

XP5L,  T51, 

451,  XWG51 

,  CF51 

» 

s 

DELTA51,  THETA51 

» 

1 

FAR51 

32 

IF  ( H 

.EQ. 

1)  IN  =  1 

WRITE 

(6, 

15  2) 

A  NUN  (IN), 

XP51,  T51 , 

H51  , 

XWG51 , 

$ 

CF51 , 

DELTA51  , 

1 

THETA51, 

FAR51 

,  ETALPT,  PR 

:.PT 

WRITE 

(6, 

153) 

A  NUN  (14), 

<  P52,  T  52 , 

H52, 

XWG51, 

$ 

CF52, 

DELTA52, 

1 

THETA52, 

FAR51 

• 

WRITE 

(6, 

154) 

ANUi  (14), 

<  ®53,  T  53 , 

H53, 

XWG53, 

2 

CF53, 

0ELTA53, 

1 

THET  A  53, 

FAR53 

WRITE 

(6, 

15  5) 

XPS3,  T53, 

4  56,  XWG6, 

CF€  , 

DELTA6 , 

$ 

THETA6 

,  FAR6, 

1 

DPE 

WRITE 

(6, 

15b) 

ALT 3 

WRITE 

(6, 

157) 

CS 

WRITE 

<6, 

158) 

DELTAH8 

WRITE 

(6, 

159) 

DEL3AHB 

WRITE 

(6, 

160) 

N 

WRITE 

(b, 

161) 

M 

WRITE 

(6, 

16  2) 

HPD3 

WRITE 

(6, 

163) 

HPEXT 

WRITE 

(6, 

164) 

HPMAH3 

WRITE 

(6, 

16  5) 

HPM 

WRITE 

(f>» 

166) 

HPA3S 

WRITE 

(6, 

167) 

HPL3C 

WRITE 

(6, 

168) 

HPHP3 

WRITE 

(b, 

169) 

XWF4 

53 


O  O 


WRiTE  (6, 

170) 

CFF 

WRITE  (6, 

171) 

HV 

WRITE  16, 

17  2) 

WLPC 

WRITE  (6, 

17  3) 

LCSl 

WRITE  (6, 

17  4) 

WHP3 

WRITE  (6, 

175) 

HCS. 

WPITE  (6, 

176) 

®CT1 

WRITE  (6, 

177) 

HPHT 

WRITE  (o, 

176) 

PCT2 

WRITE  (6, 

17  9) 

HPLT 

IF  <M  .EQ, 

,  4)  WRITE  (6, 

WRITE  (6, 

151) 

FN 

WRITE  (6 , 

182) 

FG 

WRITE  (6, 

18  3) 

FNQWA 

WRITE  (6, 

164) 

ETATHM 

WRITE  <6, 

185) 

V  6 

WRITE  (6, 

186) 

SFCJ 

IF  (IPAS  . 

NE.  2) 

GO  TO  33 

DPE 

=  0.0 

PRHPT 

=  l.'j 

PRLPT 

=  1.0 

H  V 

=  18400 

1.  0 

GO  TO  1 

T 1 

=  Ti  ♦ 

459.67 

T  2 

=  T2  * 

459.  67 

T  21 

=  T  21  ♦ 

459.57 

T  22 

*  T22  ♦ 

459.67 

TB 

=  T9  ♦ 

459.67 

T  3 

=  T3  4 

459.67 

T  4 

=  T4  4 

459.  67 

GO  TO  11 

•**  EPPOR  STATEMENTS 


34 

WRITE 
GO  TO 

<6, 

5 

187) 

35 

WRITE 

STOP 

(6, 

188) 

PPL®C, 

pfhpc,  etalpc, 

ETAH°C 

36 

WRITE 

STOP 

(6, 

18  9) 

WF40S, 

T4,  -TAB 

37 

WPITE 

STOP 

(6, 

190) 

PRHPT, 

PRLPT,  PTAHPT, 

ETALPT 

100  FORMAT  {  /  /  t  33X,  48HTHIS  PROGRAM  HAS  DEVELOPED 

$  BY  COMPONENT 

IS  BRANCH  ,  2 3H TURBINE  ENGIVE  3TVISION,  /  ,  34X, 

$  26HM0DIFIE0  BY  PO 

2WER  DIVISION,  42H  U.S.  AIR  FORCE  8ERO  PROPULSION 
l  LABORATORY  ) 
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1)1  FORMAT  <  1H1 1  2 OX,  18H**»INPUT  VALUES***  ) 

102  FORMAT  (  11A6  ) 

103  FORHAT  <  1H1  /  1H  ,  40X,  11A5  > 

104  FORMAT  (  1H0,  20X,  iaH***IN3UT  VALUES***  ) 

105  FORMAT  (  /  /  ,  5X,  lbH ALTITUDE  (FT)  =  ,  F3.0  ) 

106  FORMAT  (  5X,  17HPRE3SJRE  (PSI)  s  ,  F6.1  ) 

1)7  FORMAT  (  5X,  29H0ELT A  TEMP  -?OM  1962  ATMOS  =  ,  F8.5 
l  ) 

108  FORMAT  (  /  ,  5X,  16HBEGIN  INDICATOR  =  ,  II  ) 

109  FORMAT  (  5X»  loHENO  INDICATOR  =  ,  Ti  ) 

IIP  FORMAT  (  /  ,  5X,  2BHNO  OF  COMPRESSOR  STAGES  =  ,  11 

S  ) 

111  FORMAT  (  5X,  23HNO  OF  TURBINE  STAGES  =  ,  II  > 

112  FORMAT  (  /  ,  5X,  7HF.OM  =  ,  FID. 5  ) 

113  FORMAT  (  5X»  11HMECH  EFF  =  ,  F8.5  ) 

114  FORMAT  (  5X,  29WPA02  EXTRA0TFD  FOR  ACCESS  =  ,  F4.2 
S  ) 

115  FORMAT  (  5X,  21H0ESI GN  SHAFT  POKER  =  ,  F10.4  ) 

116  FORMAT  (  /  ,  5X,  24HIMLET  DUCT  PRESS  LOSS  =  ,  F4.2 

S  ) 

117  FORMAT  (  5X,  18HLPC  PRESS  PATIO  =  ,  F6.3  ) 

118  FORMAT  (  5X,  loHLPC  EFF  =  ,  F8. 5  ) 

119  FORMAT  (  5X,  19HlPC  SEAL  L EA <A5E  =  ,  FT.l  ) 

120  FORMAT  <  5X,  27HLPC  3LEF0,  TO  BE  USED  LATEP  ) 

121  FORMAT  (  /  ,  5X,  24HINTER-0 OMP  PRESS  LOSS  =  ,  F4.2 

$  ) 

122  FORMAT  (  5X,  24HBLEE0  DUCT  »RESS  LOSS  *  ,  F6. 3  ) 

121  FORMAT  (  /  ,  5X,  18HHPC  PRESS  RATIO  =  ,  F5.3  ) 

12«»  FORMAT  <  5X,  10HHPC  EFF  =  ,  Ff>.  5  ) 

125  FORMAT  (  5X,  19HMPC  SEAL  LEAKAGE  =  ,  F4. 2  ) 

126  FORMAT  (  /  ,  5X,  23HMAIN  COMP  PRESS  LOSS  =  ,  F8.5 

$  ) 

127  FORMAT  (  5X,  22HMAIM  COMBUSTION  EFF  =  ,  F8.5  ) 

128  FORMAT  (  5X,  27HDESIGN  FJEL  FLOW  (LB/HR)  =  ,  F6. 2  ) 

129  -ORMAT  (  5X,  45HMAIN  COMB  FUEL  HEATING  VALUE  AT  T4 
%  FOR  JP4  =  , 

1  F6.n  ) 

130  FORMAT  (  5X,  39HC0M3  DISCHARGE  TEMP  (OcGREF 
S  RANKINF)  =  ,  F7.2  ) 

131  FORMAT  (  5X,  39HHPT  INLET  NOZ  COOLTNG/HPC  INLET 
S  FLOW  =  ,  F8.5  ) 

132  FORMAT  (  /  ,  5X,  18HMPT  PRESS  RATIO  =  ,  F5.3  ) 

133  FORMAT  (  5X,  13HHPT  EFF  =  ,  r5.2  ) 

134  FORMAT  (  5X,  30HTUR3  COOLI NS7HP0  INLET  FLOW  =  , 

$  F8.5  ) 

135  FORMAT  (  /  ,  5X,  35HTNTER  TURB  DPESS  LOSS  (P5-P51) 

$  /P5  =  ,  F5.2 

1  ) 

136  FORMAT  (  /  ,  5X,  13HLPT  PRESS  RATIO  =  ,  F5.3  ) 

137  -ORMAT  (  5X,  1UHLPT  EFF  r  ,  ?e.5  ) 
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13d  FORMAT  (  5X,  39HFXHAUST  DUCT  P>*E5S  LOSS  (P*3-P1>/ 

$  P53  =  ,  F3.1  ) 

139  FORMAT  (  5X,  29HLPT  COOLING/  RPC  INLET  FLOW  =  ,  F8.5 
5  ) 

1*0  FORMAT  (  lHlf  52X,  19H***0UTPUT  VALUES***  ) 

141  FORMAT  (  1H0  ,  9X,  21HSTATION  PRtSSURF-,  A3,  2X, 

$  29HTEMP-DEG. 

IF  ENTHALPY  FLOW,  3X,  42HC0R  ;LOW  DELTA 

$  THETA  FUEL/AIR 

2  EFF,  2 X ,  18H PRESS  RATIO  OrL  p  ) 

142  FORMAT  (  1H1,  9X,  2LHSTATI ON  PRESSURF-,  A3,  2X, 

$  29HT EMP-DtG • 

IF  ENTHmL°Y  FLOW,  3X,  42HCOR  rLOW  DELTA 

$  THETA  FUEL/AIR 


2 

EFF,  2X,  18HPRFSS 

RATIO  DEL 

P  ) 

143 

FORMAT  (  1H  ,  bX, 

9HAMSI ENT 

,  3X, 

F8. 3 

S 

3 X,  F8.2, 

1 

21 X ,  F7.3,  2X,  F7.3  ) 

14*, 

FORMAT  (  1H  ,  6X , 

13HLP  COM* 

INLET 

,  4X , 

$ 

F8.2,  3X, 

1 

F8.2,  2X,  F  7 •  3 , 

3X,  F7.3,  2  X 

,  *7. 

3,  2X 

$ 

F7.5,  2X, 

2 

F7.3,  2X,  F7.5  ) 

145 

FORMAT  (  1H  ,  6X , 

12HLP  COM* 

Exrr, 

5X, 

$ 

F8.2,  3  X ,  F8.2, 

1 

2X,  F7.3,  3  X ,  F7 

.3,  2X ,  F7»  3 

,  ?x , 

F7.3 

s 

) 

14b 

FORMAT  <  1H  ,  eX, 

13HHP  COMP 

INLET 

,  4  X , 

$  F8.2,  3  X , 

1  F8.2,  2X,  F  7.  3  ,  3X,  F7.3,  2X,  -7.3,  2X,  F7.3,  12X, 
$  F  7.  5  ,  2X , 

2  F7.3  ) 


147 

FORMAT  ( 

1H  ,  bX,  A3,  1 3  ROOM  *  EXIT/,  /  1H 

,  8X, 

S 

12HBURNER 

INLET 

1 

f  3X  9  F  6  •  3 

,  5X,  F8. 2,  3 X ,  -8.2,  2X, 

F7.3,  3 X , 

F7.3, 

s 

2X,  F7 1  3 > 

?X, 

2F7.3,  2 1 X , 

F7.  3  ) 

148 

FORMAT  < 

1H  ,  oX,  11H3JRMER  EXIT, 

6X,  F8.3, 

5X, 

$ 

F8#  2  9  3Xj 

F8.  2, 

1 

2X9  F7 • 

3,  3  X ,  F7.3,  2X,  F7.3,  2X 

,  F7.3,  ?X 

,  F8.5 

S  2X ,  F7.5, 

2  11X,  F7.5  ) 

149  FORMAT  (  1H  ,  6X,  12HTU*9  1  INLET ,  5X,  F8.3,  5X, 

$  F8.2,  3X ,  F8.2, 

1  2X ,  F7.3,  3X,  F7.3,  2X,  F7,3,  2X,  F7.3,  2X,  F8.5, 
$  2X,  F7.5,  2 X, 

2  F7.3  ) 

150  FORMAT  (  1H  ,  6X,  11HTUR3  1  EXIT,  6X,  ^8.  3,  5X, 

$  F 8. 2  ,  3X,  F 8.  2, 

1  2X,  F7.3,  3  X,  F7 •  3,  2X,  F7,3,  2X,  F7.3,  2X,  F8.5, 

$  2(  X,  F 7 •  5  ) 
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151 

FORMAT  ( 

1H  , 

6X, 

11HTURB 

1 

EXIT,  / 

1H 

» 

8X, 

$ 

12HPLUS  COOLING 

> 

1 

3X  f  F8. 3 

,  5X, 

F8 

•  2» 

3X, 

FA. 

2, 

2X,  F7 

.3 

» 

3X, 

F7.3 

$ 

2 X i  F7.3, 

2X, 

2 

F7.3,  2X 

,  F8. 

5  ) 

152 

FORMAT  ( 

1H  , 

6X» 

4MTJR3, 

A3 

f 

5HINLET 

» 

5X 

»  F8 

.?» 

» 

5X,  F8.2, 

3X, 

1 

F8.2,  2X 

t  F7. 

3, 

3Xf 

r7.  3 

t  2 

•1 

^1 

• 

2X 

• 

F7.3 

,  2X 

S 

F8.5,  2X, 

2 

F7.5,  2X 

,  F7. 

3  ) 

153 

FORMAT  ( 

1H  , 

6X, 

4HTURB, 

A3 

f 

5HEXTT 

» 

5X 

»  F8 

.3, 

$ 

5X,  F8.2, 

3X, 

1 

F8.2,  2X 

,  F7. 

3, 

3X, 

F7 .  3 

,  2 

x> 

~  7 . 3 1 

2X 

» 

F7.3 

,  2X 

S 

F8.5  ) 

154 

FORMAT  ( 

1H  , 

6X, 

4MTUR3, 

A3 

9 

4MEXIT, 

/ 

1H  , 

ex. 

$ 

12HPLUS  COOLIN 

1G,  3X,  F8.3 

»  5X, 

F8 

.2, 

3X, 

F  8. 

2f 

2X,  F7 

•  3 

» 

3X  . 

F7.3 

3 

2X,  F7.3, 

2X, 

2F7.3,  2X,  F8.5  ) 

155 

FORMAT  ( 

1H  , 

6X, 

17H 

FXMT/ST 

ATIC  COND 

» 

F 

8.  3, 

5X, 

S 

F8.2,  3X, 

1 

F8.2,  2X 

.  F  7. 

3, 

3X, 

F7.  3 

,  2 

X, 

c7.3, 

2X 

* 

F7.3 

,  2X 

?  P8.5,  20X, 

2  F7. 5  ) 

15b  FORMAT  {  /  /  ,  5X,  16H ALTITUDE  (FT)  =  ,  F3. 0  ) 

157  FORMAT  (  5 X»  17HSPEED  OF  SOUMD  =  ,  F*.2  ) 

158  F ORMAT  (  5X»  29HDELTA  TEMP  FROM  1962  ATMOS  =  ,  F5.2 
$  ) 

159  FORMAT  (  5X,  30HDELTA  PRESS  FROM  1962  PRESS  =  , 

$  -4. 2  ) 

160  FORMAT  (  /  ,  5X,  23HNO  OF  0  OMP  STAGES  =  ,  II  ) 

161  FORMAT  (  5X,  20HNO  O^  TURB  STAGES  =  ,  II  ) 

162  FORMAT  (  /  ,  5X,  18HDE5 IGN  SHA-T  HP  =  ,  F10.4  ) 

163  FOPMAT  (  5X,  16HSHP  EXTRACTED  =  ,  F10. 5  ) 

164  FORMAT  (  5X,  28HSHP  cXTRACTED,  COR  FOR  SL  =  ,  FID. 3 

S  ) 

165  FORMAT  (  5X,  34HSHP  EXTRACTED.  COR  FOR  MECH  EFF  =  , 
S  F10.5  ) 

166  FORMAT  (  5X,  31HHP  EXTRACTED  FOR  ACCESSORIES  =  , 

$  F5.3  ) 

167  FORMAT  (  5X,  2QHHP  MEEDEO  FOR  L3C  =  ,  F7. 3  ) 

168  -ORMAT  (  5X,  20HHP  MEEDEO  FOR  HPC  =  ,  F15.5  ) 

169  FORMAT  (  /  ,  5X,  23HCAL DUIATED  FUEL  FLOW  =  ,  F9.4 

$  ) 

170  FORMAT  (  5X,  17HFUEL  -LOrf  COEF  =  ,  F15.*  ) 

171  FORMAT  (  5X,  45HNAIM  COM3  FUEL  MFATING  VALUE  AT  T4 

$  FOR  J34  =  , 

1  F8.1  ) 

172  -ORMAT  (  /  ,  5X,  23HWOR<  Or  L°C  ( BTU/LB)  =  ,  F6.3 

$  ) 
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173  FORMAT  <  5  X,  33HLPC  SFAL  LEAKAGE  FLOW  FRACTION  =  , 
i  F7.5  ) 

174  FORMAT  (  /  ,  5X,  23HW0R<  OF  H3C  ( BTU/LB)  =  ,  F7.  3 

$  ) 

175  FORMAT  (  bX,  33HHPC  SEA.  ..EAKASE  FLOW  FRACTION  =  , 

$  F8.5  ) 

176  FORMAT  (  /  ,  5X,  13HH?T  WORK  °ERCFNT  =  ,  F5.3  ) 

177  FORMAT  (  bX,  20HHORSEPOWER  0 r  MST  =  ,  F15.5  ) 

178  FORMAT  (  /  ,  5X,  13H.PT  WORK  PERCENT  =  ,  F5. 3  ) 

179  FORMAT  (  5X,  20HHOR5LPOWEP.  OF  L3T  =  ,  F1F.5  ) 

180  FORMAT  (  5X,  32HTHI3  IS  ASTJALLf  A  POWER  TUR9INE  ) 

181  FORMAT  (  /  ,  5X,  13HNET  THRUST  =  ,  F1Q.5  > 

182  FORMAT  <  bX,  15HGR0SS  THRUST  =  ,  F10.5  ) 

183  FORMAT  (  bX,  28HNET  THR J 5 T /E  >1 GI ME  AIRFLOW  =  ,  F10.5 

$  ) 

184  FORMAT  (  bX,  14HTHERMAL  EFF  =  ,  F8.5  ) 

185  FORMAT  (  5X,  2hHEXHAUST  EXIT  VF.OCTTY  =  ,  19X,  F9.3 

$  ) 

186  FORMAT  (  5X,  41HEN0I NE  UNINS  T  AL.ED  SFC  (LB  FJEL/HR- 
$  HP)  =  ,  F8.4 

1  ) 

187  FORMAT  (  1H0  ,  12HLP3  MOT  USE  0  ) 

lad  FORMAT  (  1H0»  44H***  ERROR  IN  COMPRESSOR  INPUT 

$  PARAMETERS  ***, 

1  /  ,  F  X |  8HPRLPC  =  ,  F7.?,  *X,  BHPRHPO  =  ,  F7.3, 

$  4  X ,  9HETALPC  - 

2,  F9.5,  4X ,  9HETAHPC  =  ,  F9,5  ) 

189  FORMAT  (  1HO ,  43H***  ER^OR  IN  COMBUSTOR  INPUT 
$  PARAMETERS  ***,  t 

1  ,  5X,  8HWF4DS  =  ,  F6,2,  ax,  bHU*  =  ,  F7.  2 ,  4X, 

S  7HETA9  =  ,  F8 .5 

2  ) 

190  FORMAT  (  1H0  ,  41H***  ERROR  IM  TURBINE  INPUT 

$  PARAMETERS  ***,  /  , 

1  5X,  8HPRHPT  =  ,  F5.3,  4  X,  7HPR'.PT=  ,  F5.3,  5X, 

$  9HETAHPT  =  , 

2  F5.2,  4X,  9HETALPT  =  ,  -5,1  ) 

END 
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SUBROUTINE  ATMOS  (HeT,  M,  DELTA) 

THIS  IS  A  SUBROUTINE  TO  COMPJTE  CERTAIN  ELFMrNTS  OF  THE 
$  1962 

U.S.  STANDARD  ATMOSPHERE  UP  TO  95  KILOMETERS. 

CALLING  SEQUENCE... 

CALL  ATMOS  (HFT,  TM,  DELTA) 

HFT  =  GEOPOTzNTIAl  A_TITJDE  (FEET) 

TM  =  MOLECULAR  SCALE  TEMPERATURE  (DEGREES 
S  RANKINE) 

DELTA  =  RATIO  OF  PRESSURE  TD  THAT  AT  SEA  LEVEL 

ALL  DATA  ANO  FUNDAMENTAL  CONSTANTS  A RF  IN  THR  METRIC 
$  SYSTEM  AS 

THESE  QUANTITIES  ARE  DEFINED  AS  EXACT  IN  THIS  SYSTEM. 

THE  RADIUS  OF  THE  EARTH  (REFT59)  IS  THE  VALUt 
$  ASSOCIATED  WITH  THE 

1959  ARDC  ATMOSPHERE  SO  THAT  PROGRAMS  CURRENTLY  USING 
$  THE  LIBRARY 

ROUTINE  WILL  NOT  REQUIRE  ALTERATION  TO  USF  THIS 
$  ROUTINE. 

D IMF NS I  ON  ALM  (10),  DELTAS  (1J>,  HB  (10),  TMB  (10) 

DATA  (H8(  I)  ,  1=1,  1C)  f  -  5.,  !).,  11.,  20., 

S  32.,  47.,  52.  , 

1  61.  ,  79.  ,  88.743  / 

DATA  (TMB (I) ,  1=1,  10)  /  320,65,  288.15,  216.65, 

$  216.65, 

1  228.65,  270.65,  27J.65,  252.65,  180.65,  180.65  / 

DATA  (DELTAB(I) ,  1=1,  10)  /  1.75363,  1., 

S  2.23361L  -  Cl, 

1  5.4C328L  -  02,  8.56663E  -  53,  1,09455-  -  03, 

$  5.822  89L  -  04, 

2  1.79718F  -  04,  1.0241s  -  05,  1.S223E  -  06  / 

DATA  (ALM  (I)  ,  1  =  1,  1C)  /  -  5.5,  -  6.5,  ,  1., 

$  2.8,  0., 

1  “  2.,  “4.,  0.,  0.  t 

DATA  RFFT59  /  2. 035553 1 Ft 7  /,  SZ  f  9.8fc6Cc  /,  AMZ 
$  /  28.9644  /, 

1  RSTAR  /  8.31432  /,  FTTOKM  f  3.048E  -94/ 

CONVERT  GEOPOTLNTIAL  ALTITUDE  TO  GEOMETRIC  ALTITUDE 
ZFT  =  HFT  *  REFT59  /  UEFT59  -  HFT) 

CONVERT  HFT  ANO  ZFT  TO  KILOMETERS 

7  =  FTTOKM  *  2=7 

H  =  FTTOKM  *  HFT 

IF  (H  .IT.  -5.0  .OR.  Z  * 51 •  90.0)  GO  TO  6 
00  1  M  =  1,  10 
IF  (H  -  HB (M) )  2,  3,  l 

1  CONTINUC 
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-•  •  t 


f  ^ 


GO  TO  6 

2  M  =  M  -  1 

3  DElH  =  H  -  HR(M) 

IF  (ALM(M)  .tQ.  C.C)  uO  TO  4 
TMK  =  TM8 <H)  *■  ALMOO  *  PE.H 

GRAOIFNl  IS  NON  ZERO,  PA  j  E  10  ,  EJ'JATTON  I.2.1G-C3) 

UElTA  =  0ELTA31H)  *  { (TM3  <M)  /  TMK)  **  (GZ  *  AMZ 

$  /  (RSTAR  * 

1  ALM(M)))) 

GO  TO  5 

4  TM<  =  TMR(M) 

GPAOItNT  IS  ZTRO,  PAGE  1J ,  EQ JATIQN  I.2.1QM4) 

DTLT  A  =  OtLTA 3 ( M)  *  EX»<  -  GZ  *  AMZ  *  QFLH  / 

$  ( PST  A  3  111 
1  TM0CM))) 

CONVERSION  TO  ENGLISH  UNITS 

5  TH  =  1.8  *  TMK 

6  RETURN 
END 
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SUBROUTINE  PROCOM  (FARX,  TEX,  CSEX,  AKFX,  CPEX, 

$  REX,  PHI,  HEX) 

CALCULATES  THERMO  GAS  PROPERTIES,  AIR  OR  JP-4/A1R 
$  MIXTURE 

COMMON  /WPAFR/  ICOJMT 
IF  <F  ARX  « LE.  C.  057623)  GO  TO  1 
FARX  =  0.067623 

WRITE  (6,  100) 

1  IF  (TEX  .GE.  3CQ •)  GO  TO  2 

TEX  =  300. 

IF  (ICOUNT  .GT.  1)  GO  TO  2 
WRITE  (6,  101) 

2  IF  ( TCX  .LE.  4000.)  30  TO  3 

TEX  =  4000. 

WRITE  (6,  102) 

3  IF  (FARX  .GE.  0.3)  GO  TO  4 

FARX  =  0.0 

WPITE  (6,  10  3) 

C  AIR  PATH 

4  CPA  =  ( ( (( ((1.3115540E  -  25  *  TEX  -  1.4526770F 

3-21)  *  TEX  I 

1  ♦  7. 6215767  £  -  18)  *  T£X  -  1. 5128259':  -  14)  *  TEX 
$  -  6. 7178376E 

2  -  12)  *  TtX  ♦  6.55194865:  -  08)  *  TEX  -  6.1536879E 

3  -  «5)  *  TEX  ♦ 

3  2. 50 20051E  -  01 

HEA  =  ((((((  (1.  264442FS  -  25  *  TFX  - 

$  2.0752522E  -  22)  *  TIX 

1  ♦  1. 270 2630E  -  18)  *  1IX  -  7.0256518E  -  15)  *  T“X 
$  -  1.  6794594F 

2  -  12)  *  TEX  +  2. 183982 3 s,  -  08)  *  TEX  -  2.5763443E 
1-05)  *  TEX  ♦ 

3  2.5)2w051L  -  01)  *  TiX  -  1.75S8886E  ♦  00 

SEA  =  +  2.5023351-  -  31  *  ALOG(TEX)  ♦ 

?  (((( ((1.4450 767E  - 

1  25  *  TEX  -  2.4211238*  -  2?)  *  rcX  ♦  1.5243153c  - 
S  18)  *  TEX  - 

2  3. 7820648E  -  15)  *  TEX  -  2,  2392793E  -  12)  *  TtX  ♦ 

S  3. 2759743E  - 

3  08)  *  TEX  -  5. 1576879F  -  05)  *  TEX  «■  4.54323J0E  - 
$  02 

IF  (FARX  • EQ«  0.0)  GO  TO  5 
3  FUEL/AIR  PATH 

CPF  =  (  (<( <(7.2b7B7l0t  -  25  *  TtX  -  1.3335668E 

$  -  20)  *  TEX 

1  ♦  1.  ( 212913E  -  16)  *  TEX  -  4.2031104E  -  13)  *  TEX 
$  ♦  9. 9666793E 

2  -  10)  *  TEX  -  1. 37719"1£  -  OF)  *  TEX  +  1.2258630E  1 

$  -  03)  *  TtX  ♦ 

i 
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3  7.  3816638L  -  02 

HEF  =  <((<((  <9.  J848388E  -  25  *  TEX  - 

3  1.9050949E  -  21)  *  TEX 

1  ♦  1.7021525E  -  17)  *  TEX  -  B.41C2208E  -  1M  *  TEX 
$  ♦  2*  ■♦921698E 

2  -  10)  *  TtX  -  4#  S9J  6332c.  -  T7>  *  TEX  6.129315PE 
$  -  04)  *  TEX  * 

3  7.  3816638E  -  i>2)  *  TEX  *  3.0581530E  ♦  Pi 

SIF  =  +  7 • 38166381  -  02  *  ALOG(TEX)  + 

$  (  (((  (  (l.C  382670t  - 

1  2E  *  TEX  -  2«  2226118E  -  21)  *  TEX  ♦  2.U425826F  - 
$  17)  *  TtX  - 

2  1.  0  512776F  -  13)  *  TEX  ♦  3.  322S928E  -  It)  *  TEX  - 
$  5«8d59505L  - 


3 

u  7)  * 

TEX 

+  1.2258630E  -  PM  * 

TEX 

+ 

6.48339*t  - 

f 

ill 

5 

CPEX 

r 

(CPA  +  FARX  *  CPF)  / 

(1. 

♦ 

FARX) 

HEX 

= 

( HFA  +  FARX  *  HE")  / 

(1. 

+ 

FARX) 

PHI 

s 

(SEA  +  FARX  *  S5F)  f 

(1. 

♦ 

FARX) 

AMW 

= 

28. 97  -  .  9461B6  *  FA*X 

REX 

= 

1 • 966375  /  AM* 

A  K£X 

= 

CPFX  /  (CPEX  -  PE  X) 

CSEX 

r 

SORT(AKEX  *  REX  •  TEX 

* 

25031.37) 

RETURN 

130 

FORMAT 

(  1H0,  b3HIN3UT  FUEL-AIR 

*ATIO 

ABOVE  LIMITS 

$ 

IT  HAS 

3EEN 

IRE  SET  TO  0.0b7623,  6HSf$S$$  ) 

101  FORMAT  (  1HU,  3FHPRJC0M  INF J T  TEMP-RATURt.  PELOW 
$  3u0« ,  6H$32$i$ 

1  ) 

102  FORMAT  (  1H0  ,  36HPR3C0M  INPUT  TEMPERAT'JRt  AROVF 
$  4C0C. i  6HSSS$3S 

1  ) 

103  FORMAT  (  1H0  ,  36HPR0C0M  INPJT  FJFL-AIR  RATIO  9F.L0W 
S  ZERO, 

1  6Hl$S$S$  ) 


SUBROUTINE  THERMO  (*X,  MX,  TV,  3X,  AMX,  CPT,  GA,  L, 
S  FAR,  K> 

C  ***  THEPMODUNAMIC  PROPf-kTIES  BASED  ON  PROCOH 
COMMON  /WPAFB/  ICOUNT 
ICOUNT  =  1 

FX  =  A. 

IF  (L  .EO.  1>  FX  =  FAR 

IF  <K  .EQ.  II  GO  TO  t 

CALL  PkOCOM  (FX,  TX,  CS,  AK,  CP,  R ,  PHI,  HX) 

GO  TO  3 

1  TX  =  U,  *  HX 

00  2  I  =  1,  15 

CALL  PROCOH  (FX,  TX,  CS,  A<,  C*,  P,  PHI,  Ml 

ICOUNT  =  ICOUNT  ♦  1 

DELH  =  HX  -  M 

IF  (ABS  (DELH)  .LE.  O.OOCOIMX)  SO  TO  3 

2  TX  s  TX  ♦  4*  *  DFLH 

WRITE  (6,  100) 

3  SX  *  °HI  -  R  *  ALOS(PX) 

AMX  *  1.936375  /  R 

CPT  =  CP 

GA  =  AK 

RETURN 
C 

100  FORMAT  (  31HONO  CONVERGENCE  IN  THTOHOtS$t$S  ) 

3 
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SUBROUTINE.  COMo  (PIN,  TOJT,  HIN,  FTA3,  0P8,  POUT» 

$  HOUT,  SOUT, 

1  FAROUT,  HV) 

t*»  COMPUTES  STOICHIOMETRIC  T t  NPER  AT'JRF  IN  COMBUSTOR  -  If 
$  HEATING 

VALUE  FOR  FUEL  ANO  W"4DS  4RE  NOT  INPUT*  T4,  FTA8* 

$  OPB  ARE  INPUT, 

DIMENSION  X  (9) 

X  ( 2)  =  0,0 

X  ( 3)  =  0.3 

POUT  =  PIN  *  (1,3  -  CPP) 

HV  =  ((((((  -  ,  4594317  £  -  19  *  TOUT)  - 

$  .20341 ItoF  -  15)  * 

1  TOUT  «■  ,27ft  3C43  c  -  11)  *  TDJT  f  .  2051501E  -  07)  * 

$  TOUT  - 

2  •  245 3116L  -  l3)  *  TOUT  -  .9433296'";  -  01)  *  TOUT  + 
S  .  1845537F  + 

3  05 

CALl  THERMO  (POUT,  HOIJTA,  TOJT,  XI,  X2,  X3,  X4,  0, 

S  0.0,  C) 

FAROUT  =  ( HOUTA  -  HIN)  (  ( (HV  *  ETA8)  -  (HOUTA  - 
$  HIN) ) 

FARS  =  FAROUT 

CALL  THERMO  (POUT,  HDUT ,  TOUT,  3DJT,  Tl,  Y2,  Y3,  l, 

$  FAROUT,  C) 

1  DFLFAR  =  (FARS  -  FAROUT  >  f  FAROUT 

CALL  A F QUIP  (X,  TOUT,  DtLrAR»  0.,  15.,  .0201,  .95, 

$  TNEW,  ICON) 

IF  (ICON  -  2)  2,  4,  3 

2  TOUT  =  TNFW 

HV  =  ((((((  -  .45943171  -  19  *  TOUT)  - 

$  .203tll6F  -  15)  * 

1  TOUT  +  ,278  3643;-  -  11)  *  TOUT  f  .2V51501E  -  07)  * 
$  TOUT  - 

2  . 245311DE  -  03)  *  TOUT  -  .94332961  -  ID  *  TOUT  ♦ 
$  , 18455 37F  + 

3  A5 

CALL  THERMO  (POUT,  HDUTA,  TOUT,  XI,  X2,  X3,  Xh,  0, 

$  G  •  C  $  C  ) 

FAROUT  =  (HOUTA  -  HIN)  /  ( (HV  *  CTAB>  -  (HOUTA  - 
$  HIN) ) 

GO  TO  1 

3  WRITE  (6,  10  0) 

4  CALL  THERMO  (POUT,  HOUT,  TOUT ,  SOUT,  Xi,  X2,  X3,  1, 
$  FAROUT,  0) 

RETURN 


IjQ  FORMAT  (  46H  STOICHIOMETRIC  TEN*  WTLL  NOT  CONVERSE 
$  IN  COMB  ) 


64 


SUBROUTINE  C OH 32  (PIN,  TDUT ,  HI*,  ’•TAB,  DPB,  POUT, 

I  HOUT,  SOUT, 

1  FAR4,  HV) 

**♦  SIMILAR  TO  COMB  -  HEATING  VALJE  FDR.  pUEL  IS  OPTIONAL, 
USSR  MUST  INPUT  T<*,  iTAB,  DPB,  Ws4DS. 

DIMENSION  X  (9) 

X (2)  =  0.0 

X  (3)  =  0.0 

POUT  *  PIN  *  <1.3  -  DPP) 

1  CALL  THERMO  (POUT,  H3UTA,  TOJT,  XI,  X2,  X3 ,  X4,  0, 

$  0.0,  0) 

FARCUT  =  (HOUTA  -  HIN)  t  ((Hi/  *  tTAB)  -  (HOUTA  - 
J  HIN)  ) 

FARS  =  FAROUT 

DFLrAR  =  ( F AP,4  -  ?ARS)  (  r  AR4 

DIR  *  SORT (FAR*  /  FSRS) 

CALL  AFQUIP  (X,  TDUT,  DEL-AR,  0.,  20.,  .0^01,  DIR, 

1  TNEW,  ICON) 

IF  (ICON  -  2)  2,  4,  3 

2  TOUT  =  TNEW 

GO  TO  i 

3  MPITL  (6,  100) 

4  CALL  THERMO  (POUT,  HDUT ,  TOUT,  SDUT,  XI,  X2,  X3,  1, 

$  F AKA  ,  G) 

RrTURN 

loO  FORMAT  (  47H  STOICHIOMETRIC  TtM®  WILL  NOT  CONVERGE 

2  IN  COMB2  ) 


O  O  O  O  c : 


SUBROUTINE.  AFQUIR  (X,  AIND,  3E»END,  ANS,  A J ,  TOL, 

$  DIR,  ANEW, 
i  ICON) 

♦**  QUADRATIC  CONV ERGFNCE  ROUTINE.  ^DR  INTERPOLATION 

X ( 1) =NAME  OF  ARRAY  TO  USE 

AIND=INOt PENDANT  VARIABLE 

Oc.PENO=  DEPENDANT  VARIABLE 

ANS= ANSWE R  UPON  WHICH  TO  CONVERGE 

AJ=MAX  NUMBER  OF  TRYS 

T  OL=PERGE  NT  TOLERANCE  FOR  CONVZRGENCE 

DIR=DIF  ECTI ON  ANO  PERCENTAGE  FDR  FIRST  GUESS 

ANEW=CALCUL ATEO  VALUE  OF  NEXT  TRY  AT  T NQrPENOANT  VARIABLE 


ICON=CONTROL 


X(2)=COUNTER 
X  { 3 ) =CHOOS£S  METHOD  OF 
X ( 4) =THIRD  DEPEND  VAR 
X ( 5 ) =TH1k 0  I ND  VAR 
X(6)=SrC0ND  DEPEND  VAR 
X  ( 7 ) =SECONO  IND  VAR 
X ( 8) =FIRST  DEPEND  VAR 
X  ( 9) =FIRST  IND  VAR 
X ( 3 )  MUST  BE  ZERO  UPON 


=  1  GO  THRU  _ ODP  AGAIN 
=  2  YOU  HAVE  REACHED  ME  ANSWER 
=3  COUNTER  HAS  HIT  LIMIT? 
STORAGE 

CONVERGENCE 


FIRST  ENTRY  TO  ROUTINE 


DIMENSION 

X  (9) 

Y  = 

0, 

IF  (ANS  .uQ. 

0)  GO  TD  1 

Dfc-P 

DEPEND  -  ANS 

TOLANS  = 

GO  TO  2 

TOL  *  ANS 

1 

OEP 

DEPENO 

TOLANS  = 

TOL 

2 

IF  (ABS(OEP) 

•LE.  TO. ANS) 

IF  ( X  ( 2)  - 

A  J)  b  ,  Of  5 

3 

ANEW  = 

AINO 

X  (2) 

0. 

ICON  = 

RETURN 

2 

4 

ANEW  = 

Y 

X  (2)  = 

X  ( 2)  «■  1. 

ICON  = 

RETURN 

1 

5 

ANEW  = 

Y 

X  (2)  = 

0. 

ICON 

RETURN 

3 

6 

IF  ( X ( 3 )  .GT.  0)  GO  TO  l" 

GO  TO  3 
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7 

X  (3) 

= 

i. 

XC8) 

DEP 

X  (9) 

* 

AIND 

IF  (AIND) 

8 1  9.  8 

8 

Y 

GO  TO  4 

S 

DIR  *  AIND 

9 

Y 

GO  TO  4 

D  IP 

10 

IF  (X (3) 

- 

1.)  11.  lit 

«** 

LINEAR  GUFSS 

11 

X  (3) 

r 

2. 

X  (6) 

r 

DEP 

X  (7) 

x 

AINO 

IF  f. X ( 8 )  .EQ.  X(6)  .OR.  X  ( 9 )  .EQ.  X(7>>  GO  TO  7 
A  =  ( X(9)  -  X  ( 7)  >  /  (  X(3)  -  X(6>> 

r  =  x (9)  -  a  *  xm 

IF  (ABS(1G.  *  X (9)  )  -  A3S(Y)>  7,  7,  4 

C  ***  QUADRATIC  GUFSS 
12  X (4)  =  OEP 

X  (5)  =  AIND 


IF 

(  X  ( 7 ) 

.  EQ. 

X  ( 5) 

•  AND. 

X  ( 6) 

•  EH. 

X (4)  ) 

GO 

TO 

7 

IF 

( X  (  7  ) 

•  NE. 

X  (5) 

•  AND. 

X  ( 6) 

•  NE. 

X  (4)  ) 

GO 

TO 

13 

GO 

TO  11 

IF 

( X  ( 9) 

•  uQ. 

X  (5) 

•  AND. 

X  ( 8) 

•  E3. 

X(4)  ) 

GO 

TO 

14 

IF 

( X  ( 9) 

•  NE. 

X  (5) 

•  AND. 

X  (  8) 

•  NF. 

X  (4)  ) 

GO 

TO 

15 

X  (91  =  X  (7) 

X  (8)  =  X  (6) 


GO  TO  11 

14  X  (9)  =  X (7) 

X ( 8)  =  X (6) 

X (3)  =  1. 

IF  (X  ( 9)  )  8,  9,  8 

15  F  =  (  X(6)  -  X  ( 4)  )  /  CX(7)  -  XC5)) 

A  =  ( X  (8)  -  X  ( 4)  -  F  *  ( X  (  9)  -  X  (  F  )  )  )  / 

$  < (X ( 9)  -  X (7)  )  * 

1  (X  (  9)  -  X  (5 ) ) ) 


B 

=  F  -  A  *  ( X (5) 

«■ 

X  ( 7)  ) 

C 

*  X (4)  ♦  X(5)  * 

(A 

*  X(7) 

- 

IF 

(A 

.  EQ>  0  .AND.  B  .  E]i 

0) 

30  TO 

5 

IF 

(A 

•NE.  0  .AND.  B  .ME. 

0) 

30  TO 

16 

IF 

(A 

•  NE.  C  .AND.  3  .NE. 

0) 

GO  TO 

18 

Y 

=  -  C  /  B 

GO 

TO 

27 

IF 

(C 

•  NE.  0.0)  GO  TO  17 

Y 

=  0, 

GO 

TO 

27 

G 

=  -  C  /  A 

IF 

(G 

•IE.  O.J)  GO  TO  5 

Y 

=  SQRT(G) 

YY 

=  -  SQRT(G) 
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GO 

TO  22 

18 

IF 

(C  .NE.  0.0  GO  TO  19 

Y 

=  -  B  /  A 

YY 

=  0. 

GO 

TO  22 

19 

0 

=4.*A*C/3** 

? 

IF 

(1.  -  D)  11,  20,  21 

2u 

Y 

=  -  9  /  (2.  •  5) 

GO 

TO  27 

21 

r* 

=  SORT  (1.  -  3) 

Y 

=  (  -  B  ✓  (2.  *  A)) 

<1.  ♦ 

E) 

YY 

=  (  -  B  /  (2.  *  A>> 

» 

(1.  - 

t) 

22 

J 

=  4 

00  23  I  =  6,  8,  ? 

IF  (ABS(X(4))  .LE.  A  33  C  K  ( I> )  > 

SO  TO 

23 

23 


J  =  I 

CONTINUE 

<  =  J  +  1 

IF  <  (X  (K)-Y)MX(K)-YY)  .LE.  3.")  GO  TO  24 
IF  (ABS(X(K)-Y)  .LE*  A3S(X(K>-YY>>  GO  TO  27 

Y  =  YY 

GO  TO  27 

IF  (J  .GE.  6)  GO  TO  25 
JJ  =  J  ♦  2 

KK  =  <  ♦  2 

GO  TO  26 

JJ  =  J  -  2 

KK  =  K  **  2 

26  SLOPE  =  ( X  (KK)  -  X<<>>  /  (XUJ)  -  X(J>) 

IF  ( SL  OPt*X  <J)*(X(K)"Y)  •  ST  •  3.1)  GO  TO  27 

Y  =  YY 

X  ( 9)  =  X  (7) 

X  (8)  =  X  (6) 

X (7)  =  X (5) 

X (6)  =  X (4) 

GO  TO  4 
ENO 


24 


25 


27 
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